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ABSTRACT 

We present the rest-frame optical luminosity function (LF) of red sequence galaxies in 16 clusters at 
0.4 < z < 0.8 drawn from the ESO Distant Cluster Survey (EDisCS). We compare our clusters to an 
analogous sample from the Sloan Digital Sky Survey (SDSS) and match the EDisCS clusters to their 
most likely descendants. We measure all LFs down to M ~ M* + (2.5 — 3.5). At z < 0.8, the bright 
end of the LF is consistent with passive evolution but there is a significant build-up of the faint end 
of the red sequence towards lower redshift. There is a weak dependence of the LF on cluster velocity 
dispersion for EDisCS but no such dependence for the SDSS clusters. We find tentative evidence 
that red sequence galaxies brighter than a threshold magnitude are already in place, and that this 
threshold evolves to fainter magnitudes toward lower redshifts. We compare the EDisCS LFs with the 
LF of co-eval red sequence galaxies in the field and find that the bright end of the LFs agree. However, 
relative to the number of bright red galaxies, the field has more faint red galaxies than clusters at 
0.6 < z < 0.8 but fewer at 0.4 < z < 0.6, implying differential evolution. We compare the total light in 
the EDisCS cluster red sequences to the total red sequence light in our SDSS cluster sample. Clusters 
at 0.4 < z < 0.8 must increase their luminosity on the red sequence (and therefore stellar mass in red 
galaxies) by a factor of 1 — 3 by z = 0. The necessary processes that add mass to the red sequence in 
clusters predict local clusters that are over-luminous as compared to those observed in the SDSS. The 
predicted cluster luminosities can be reconciled with observed local cluster luminosities by combining 
multiple previously known effects. 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: clusters: general — galaxies: 
luminosity function, mass function 
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1. INTRODUCTION 

Most of the stellar mass in the local universe is 
contained in "red and dead" galaxies, i.e. galaxies 
which have stopped forming stars at an appreciably 
level and whose light is thus dominated by old, red 
stars(Hog g et al.ll2002[ ). To understand how stars form 
and galaxies are assembled, we therefore need to de- 
termine how the red galaxy population evolves through 
time. Red galaxies are located on a tight sequence in 
color and magnitude, the "red sequence" (e.g. de Vau- 
couleurs 1961; Visvanathan & Sandage 1977), and the 
very small intrinsic scatter in color implies that the red 
colors result from uniformly old stellar ages (e.g. Bower, 
Kodama & Terlevich 1998). Old ages for red sequence 
galaxies are also found by studies of their stellar indices 
(e.g. Trager et al. 1998). Some studies even find a stel- 
lar mass dependence in the mean stellar age, such that 
lower mass galaxies formed their stars at later epochs 
than those that are more massive (e.g. T homas et al. 
2005) , but this result is still controversial as lTrager et al.l 
()2008l ) find no such trend in their studies of Coma cluster 
early types. 

At face value direct lookback observations may sup- 
port these local archaeological studies as the total stellar 
mass on the red sequence may have doubled since 2^1 
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(iBell et alJ[200l iFaber et al.l[2007t iBrown et al.l [200l. 



iCimatti et alJ (|2006l ) and lBrown et alJ (|2007f ) concluded 
that this mass growth comes primarily from the addi- 
tion of low mass galaxies to the red sequence at late 
times, with the most luminous red sequence galaxies 
(L > 4L*) appearing to have been in place since z > 1. 
As iTrager et al.l (|2008D point out, however, it may be 
hard to relate the direct lookback results to studies of 
local galaxies, as the latter may be susceptible to very 
small amounts (a few percent) of late star formation. 

It is impossible to study the evolution of red galaxies 
without examining t he influence of environment . Going 
all the way back to iHubble &: Humasonl ()1931h it has 
been known that there are significant correlations be- 
tween color and environment, star formation rate (SFR) 
and star formation history (SFH) and environment, and 
morphology and environment (e.g. Dressier 1980), such 
that dense environments, e.g. the centers of galaxy clus- 
ters, have much higher fractions of red sequence galaxies 
than the field. Local studies suggest that luminous el- 
lipticals in galaxy clusters have older s tellar ages than 
those in the field (|Thomas et al.l [20051 ) but studies at 
high redshift detected no differ ence in the ages of field 
and cluster elliptical galaxies (Ivan der Wei et aP [20051 : 
Ivan Dokkum fc van der Marell l2007f ). Nonetheless, the 
large differences between clusters and the field even at 
intermediate redshift, which are measured in terms of 
the morphological fraction (e.g. Postman et al. 2005; 
Smith et al. 2005; Desai et al. 2007) and the fraction of 
star forming galaxies (e.g. Poggianti et al. 2006) implies 
that a galaxy's evolutionary path might be strongly af- 
fected by the environ ment in which it li ves a s it evolves 
through cosmic time. iPoggianti et al.l (j2006f ) postulate 
that massive ellipticals in clusters may have been formed 
at high redshift but that lower luminosity red galaxies 
are added to the cluster at z < 1. 

From a theoretical standpoint, some models (e.g. De 
Lucia et al. 2006) predict that stars in red galaxies were 
formed at high redshift and that the formation epoch of 
the stars is earlier for higher mass galaxies. It is nonethe- 
less not clear if these models can be reconciled in detail 
with the observed evolution in the increase of mass on 
the red sequence at z < 1. Also not clear is if the prop- 
erties of galaxies as a function of environment are being 
properly treated in some models as none of the commonly 
implemented processes, e.g. ram-pressure stripping, ha- 
rassment, strangulation, can reproduce the observed de- 
pendence of the red and blue galaxy fraction on e.g. halo 
mass and central halo galaxy type at low redshift (e.g. 
Weinmann et al. 2006) or at z - 1 (jCoil et al.ll2008f ). 

One way to study the evolving galaxy population is to 
use the luminosity function (LF; see Binggeli, Sandage, & 
Tammann 1988 for a review), which describes the num- 
ber of galaxies per unit luminosity. The LF encodes in- 
formation about the efficiency of star formation and feed- 
back in galaxies and how galaxies populate their parent 
dark matter halos. 

Enabled by large surveys at low redshift such as 2dF 
Folkes et al.l [l9 99l and the Sloan Digital Sky Survey 



York et al.ll2000f ) it is now possible to construct the de- 



tailed LF of low redshift galaxies in a range of environ- 
ments. For example, using the 2dFGRS. [de Propris et all 
|2003i ) measured the composite LF in a set of local galaxy 
clusters and found that clusters have a brighter charac- 



teristic luminosity and steeper faint-end slope than the 
field, with the largest difference being found for spec- 
troscopically identified non-starforming galaxies. The 
availability of these well-characterized local LF determi- 
nations provide well established reference points against 
which to measure evolution in the cluster galaxy pop- 
ulation. Simultaneously, the recent availability of deep 
multi-color photometry of intermediate and high redshift 
clusters with extensive spectroscopic follow-up have al- 
lowed the galaxy population to be studied out to z ^ 1 
in the Universe's densest regions. 

De Lucia et al. (2004; hereafter DL04) were the first to 
measure the evolution of the red sequence LF in clusters 
at high redshift by studying the ratio of luminous-to- 
faint red sequence galaxies Mum/A'faint in four clusters 
at z ~ 0.7 drawn from the ESO Distant Cluster Survey 
(EDisCS). They found that this ratio was significantly 
higher in the high redshift clusters than in the Coma clus- 
ter. Subsequ ently, this redshift tren d in Mum/A^fa i nt was 
confirmed bv lGoto et all (j2005') and' Tanaka et all (j2005l ) 
in a fe w clusters, and by D e Luc ia et al. (2007 ; hereafter 
DL07). [StotFetan ([200l . and iGilbank et al.l (^ 8) in 
sign ificantly larger sampl es. iTanaka et al.l ( 20051 ). DL07. 
and IGilbank et all (j2008D also found that the evolution 
of A^ium/A^faint depends weakly on cluster velocity disper- 
sion and DL07 and Gilbank et al. (2008) find that poorer 
systems evolve marginally slower than richer systems at 
0.4 < z < 0.8. The behavior in iTanaka et al.( (|2005[) is 
based on only one cluster and is harder to generalize. 
Tracing the evolution to z = 0, however , there i s sorn e 
disagreement between DL07 and Gilban k" et al.l (j2008f ). 
In DL07 it appears that the low-dispersion systems have 
converged to the Mum/A^faint value of the Coma cluster 
while the high-dispersion systems require significant evo- 
lution to reach the value from SD SS or Coma. On the 
other hand, the poor systems of IGilbank et al.l (j2008l ) 
have systematically higher Aium/A'faint values than rich 
systems at 0.4 < z < 0.6 and therefore need to evolve 
more at z < 0.4 to come into agreement with the local 
value. The origin of this apparent discr epancy is hard 
to track down since DL07 and Gilbank et^ ([200l use 
different effective velocity dispersion cuts and different 
magnitude limits defining the split between faint and lu- 
minous galaxies. At the same time (|Andreonll20n6L[200l 
claim little trend in A^ium/Afaint with redshift and no 
trend with velocity dispersion. In their Figure 4, how- 
ever, the amount of redshift evolution appears similar 
to that from DL07. It is also not easy to compare the 
trends w ith velocity dis persion between the two works 
since the lAndreonI (|2008f ) sample contains no clusters be- 
low 600 k m/s, which cq r nprise s a large fraction of the 
DL07 and IGilbank et all (|2008f) samples. 

This paper makes a series of advances over previ- 
ous works by computing the full red-sequence LFs from 
EDisCS and comparing them to the local red sequence 
cluster LF as determined from the SDSS. The EDisCS 
sample is the largest sample that probes well past z = 
0.5, all the way out to z = 0.8, has deep multi-band pho- 
tometry that enables the construction of rest-frame opti- 
cal luminosity functions, and has a large range in cluster 
velocity dispersion. In this paper we extend the work 
of DL07 significantly by measuring the non-parametric 
LF, fitting Schechter functions, and measuring the de- 
tailed evolution of red sequence galaxies. In doing so we 
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pay specific attention to the ability to determine mem- 
bership from galaxies with only photometry. Our large 
range in velocity dispersion permits us to study how evo- 
lution in the LF depends on velocity dispersion and our 
deep photometry make us complete well below M*. We 
also make the first comparison of the composite cluster 
red sequence LF to that in the field and measure their 
comparative evolution. This test is crucial as it spans the 
full range of galaxy environment and speaks directly as 
to whether the cluster and field red galaxy populations 
are built up at different rates. Finally we measure the 
evolution of the total light on the red sequence in clus- 
ters and discuss its implications for how mass is added 
to the cluster red sequence over time. We do not ad- 
dress in detail the total LF or that of blue galaxies as we 
show in §4.51 that LFs from photometric data can only 
be robustly computed for red galaxies. 

In this paper we examine the rest-frame optical LF of 
the red galaxies in EDisCS clusters. The rest-frame near 
infrared LF and stellar mass function will be presented in 
Aragon-Salamanca et al. (in preparation). In fJ2]we dis- 
cuss the survey strategy and describe the data. In ij3]we 
discuss our techniques for determining cluster member- 
ship. In 12] we describe our estimation of rest-frame lu- 
minosities and present our construction of the rest-frame 
optical LF. We present our results in §5, discuss them 
in §6, and summarize and conclude in §7. Throughout 
we assume "concordance" A-dominated cosmology with 
flu = 0.3, 17a = 0.7, and Ho = 70 hyo km s'^ Mpc'^ 
unless explicitly stated otherwise. All magnitudes are 
quoted in the AB system. 

2. OBSERVATIONS AND DATA 

2.1. Observations and Survey Description 

The survey strategy and description are presented in 
detail in White et al. (2005; hereafter W05) who also 
present the optical photometry and the construction of 
photometric catalogs. The near Infrared (NIR) photom- 
etry will be presented in (Aragon-Salamanca et al. in 
preparation). The spectroscopic data are presented in 
[Halliday et al. (200 4) for the first five clu s ters w ith full 
spectroscopy and in iMilvang- Jensen et al.l ()2008D for the 
full EDisCS sample. The survey description and data 
will be summarized briefly below. 

The original goal of EDisCS was to study in detail a set 
of 10 clusters at z 0.5 and 10 at z ~ 0.8. Our survey 
draws on the optically selected sample of clusters from 
the LCDCS (jGonzalez. Zaritsky. Dalcanton. fc NelsonI 
1200 It) . After confirming the presence of a galaxy surface 
overdensity at the expected position and the presence of 
a red sequence using short images with the F0RS2 in- 
strument on the VLT, we initiated deep imaging of 10 
clusters in each redshift bin. We observed every field in 
either the B, V, I, and Xg-bands or in the V, R, I, J, 
and if^-bands depending on whether the LCDCS redshift 
estimate of the cluster was at 0.5 or 0.8 respectively. The 
optical data were all obtained with F0RS2/VLT and the 
NIR data were obtained with the SOFI instrument on the 
NTT. 

From the first reduction of our imaging data we com- 
puted photometric redshift s to get a more pr ecise redshift 
estimate for the clusters ( Pello et al.ll2009h . These red- 
shifts were used to target objects for spectroscopic obser- 
vations with F0RS2/VLT. Now complete, our extensive 



spectroscopic observations consist of high signal-to-noise 
(S/N) data for ^ 30 — 50 members per cluster and a 
comparable number of field galaxies in each field down 
to / ~ 22. As explained in W05, deep spectroscopy was 
not obtained for two of the EDisCS fields (CL1122.9- 
1136 and CL1238. 5-1144), the former of which showed 
no evidence for a cluster in the initial, shallow spectro- 
scopic observations. These clusters will not be used in 
this study, leaving 18 of which one (CL1119. 3-1129) does 
not have any NIR data. 

2.2. Catalog Construction and Total Flux Measurements 

We measured two types of magnitudes for our galax- 
ies, matched aperture magnitudes and SExtractor AUTO 
magnitudes. The former are used for measuring colors 
and the spectral energy distributions (SEDs) used to fit 
the photometric redshifts. The latter are used to esti- 
mate the total magnitude of the galaxies in question. 
We describe each in turn. All mag nitudes have been cor - 
rected for galactic extinction from iSchlegel et all (fl998^. 

Before the measurement of matched aperture fiuxes, 
all images with seeing better than FWHM=0'.'8 were 
convolved to FWIIM=0'.'8. The seeing across all bands 
ranged from 0'.'6 to I'.'O with most observations having 
FWHM< 0'.'8. 

Flu x catalogs were created u sing the SExtractor soft- 
ware (jBertin fc ArnoutsI Il996h in the two image mode, 
detecting in the unconvolved (i.e. natural seeing) /-band 
image and measuring fiuxes in matching apertures in all 
other bands. Colors were measured with the same aper- 
ture in all bands, using either isophotal apertures defined 
from the detection images for those galaxies that were 
not crowded or using circular apertures with r = I'.'O 
for those galaxies that were crowded. With this dual 
choice of matched apertures we obtained a high S /N mea- 
surement of the color while minimizing the biases due to 
neighboring objects. 

Obtaining accurate total magnitudes is important 
when characterizing the LF. A true total magnitude es- 
timate is not possible, however, due to uncertainties in 
the galaxy profile at large radii coupled with an uncertain 
knowledge of the sky level. As described in W05, there- 
fore, we attempted to measure pseudo-total magnitudes 
(called "total" magnitudes throughout) in the /-band us- 
ing the AUTO magnitude from SExtractor. These mag- 
nitudes were measured on the images at their natural 
seeing. The SExtractor AUTO measurement is executed 
within a Kron like aperture (Kron 1980) and measures 
the flux within a radius corresponding to 2 times the 
first moment of the light distribution. The AUTO mag- 
nitudes for each object have a minimum aperture radius 
of 3.5 pixels (or 0'.'7). The AUTO aperture is quite large 
for bright objects but for faint objects the AUTO aper- 
ture shrinks its size to the minimum allowable limit. In 
this regime, light will be lost out of the aperture even 
for point sources, since the stellar PSF throws significant 
amounts of light beyond this minimum aperture. Such 
an effect was also noted in the absolute magnitude esti- 
mates of lLabbe et al.l (|2003l ) and we adopt their approach 
for correcting for this effect, which we summarize here. 
Correcting for this offset explicitly is difficult because 
we don't know the intrinsic profile of the galaxies whose 
photometry we wish to measure. However, a conservative 
and necessary correction can be made by accounting for 
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the light that would be missed assuming that the object 
is a point source. While the amount of light lost may be 
larger for extended objects, this robust correction must 
be made regardless of the intrinsic object shape. Since 
we only define the total magnitude consistently from the 
/-band image, and use this to scale our rest-frame lu- 
minosities (as measured in the matched apertures) to 
total luminosities, we only calculated the aperture cor- 
rection for the /-band image. This neglects the effects of 
large color gradients, but the resultant error in the total 
magnitudes should not dominate our uncertainties. We 
determined an empirical stellar curve of growth for each 
image using a set of bright, unsaturated, and isolated 
stars. Using the curve of growth, we computed the cor- 
rection as a function of AUTO aperture area and apply it 
to the AUTO magnitudes. The corrected magnitudes be- 
come our "total" magnitudes, /tot- For the two clusters 
with the worst and best seeing in the /-band we plot the 
dependence of these corrections and the AUTO aperture 
size on the /tot in Figure [1] The corrections range from 
median values of ~ 0.04 magnitudes at 20.4 < /tot < 22.4 
to 0.09 magnitudes at 24.4 < /tot < 24.9. 

To check how well this aperture correction does in 
retrieving the true total magnitude, we compared the 
/tot values to those derived from 2D profile fits to the 
/-band data using the GIM2D software (Simard et al. 
2002; Simard et al. in prep). We fit bulge-t-disk mod- 
els to the galaxies and extrapolated the profiles to get 
total magnitude estimates /gim2D- For sources with 
no nearby neighbors, /gim2D should be relatively free 
of bias (Haufilcr et al. 2007). At 20.4 < /tot < 22.4 
and 22.4 < /tot < 24.4 we find a median difference 
-^tot — ^GiM2D of 0.02-0.04 and 0.06-0.1 respectively, such 
that_/GiM2D i s sys tematically brighter. However, in 
I Simard et al.l lj2002f ) those authors used extensive sim- 
ulations to show that the GIM2D magnitudes are biased 
brighter by the same order as our measured difference 
between /tot and /gim2D, implying that our /tot magni- 
tudes indeed are good approximations to the true total 
magnitude. 

We have verified that our results do not depend sen- 
sitively on the value of the correction, as it only signif- 
icantly effects the very faintest galaxies in the sample, 
which do not dominate any of the observed effects. 

3. DETERMINING CLUSTER MEMBERSHIP 

At intermediate redshift, the contrast of a cluster 
against the background and foreground is very low and 
an estimation of the cluster galaxy LF necessitates that a 
sample of cluster members be assembled which has been 
cleaned of foreground and background interlopers. Spec- 
troscopy is obviously the most accurate method of ac- 
complishing this and spectroscopic redshifts can be de- 
termined for single objects down to / '--^ 24.5 with the 
use of 8-meter class telescopes. Nonetheless, determining 
redshifts for large numbers of cluster members in multi- 
ple clusters, even with large time allocations on 8-meter 
class telescopes, is limi ted to relatively bright magni- 
tudes, e.g. / < 22 - 23 (|Tran et al.l 120071 : iHallidav et aH 
[200l . To deter mine cluster membership for magnitude 
selected samples down to / ^ 25 it is therefore nec- 
essary to use alternate techniques. We have employed 
two methods to accomplish this, one based on photo- 
metric redshifts Zphot and one based on statistical back- 
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Fig. 1. — An illustration of how our aperture correction depends 
on apparent magnitude for two clusters with the worst and best 
7-band image quality in our sample, CL1018. 8-1211 and CL1054.7- 
1245 respectively. The x-axis in all plot is our "total" estimate of 
the /-band magnitude, hot, which is the AUTO /-band magnitude 
with a point source aperture correction. The bottom row of panels 
shows how the correction depends on hot and the top panel shows 
how the circularized AUTO aperture radius depends on hot- Only 
objects with no evidence of crowding have been used. Stars are 
indicated by blue stars. At every magnitude, the objects with 
the smallest apertures receive the largest correction. The smallest 
apertures correspond to those for stars. 

ground subtraction. These membership techniques have 
also been used in previous works on the EDisCS clusters, 
e.g. DL07, DL07. LFs computed from photometric red- 
shifts and statistical background subtraction will here- 
after be referred to as LFzp and LFss respectively. We 
discuss these two methods in this section, while in i i4.5l 
we compare LF^p and LF^g to determine the robustness 
of our results. 

3.1. Photometric Redshifts 

In general, photometric redshift techniques estimate 
the redshift of a galaxy by modeling the broad-band SED 
with a set of template spectra (e.g. Fernandez-Soto, 
Lanzetta, and Yahil 1999; Rudnick et al. 2001). The 
resulting of the template fit as a function of redshift 
gives an estimate of the redshift probability distribution 
P{z) and hence the most likely redshift. As an example 
applicat ion of phot o metri c redshift techniques to cluster 
studies, iToft et all (|2004D used their Zphot estimates to 
determine membership by taking a very wide Az = ±0.3 
slice in redshift and selected every galaxy within this slice 
as being a cluster member. A slice this wide, however, is 
~ 100 times larger in velocity than the expected veloc- 
ity width of the cluster, implying a large contamination 
from field galaxies. Also, the performance of photomet- 
ric redshifts is expected to depend on galaxy SED shape, 
e.g. blue star- forming galaxies have weak Balmer/4000A 
breaks which result in weaker photometric redshift con- 
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straints and possible larger systematic errors. This color 
dependence on the Zphot accuracy can only be quantified 
by using a large number of spectroscopic redshifts that 
span a large range of SED shape/color in the desired red- 
shift range, preferably with identical photometry. Until 
now, such large spectroscopic samples in intermediate 
redshift cluster fields have not been available. 

We explore an alternative photometric-redshift-based 
interloper subtraction technique with EDisCS, which 
tries to mitigate the disadvantages mentioned above. 
The photometric redshifts for the EDisCS sample, their 
performance, and th eir use to isolate c luster members, is 
described in detail in lPello et al.l (|2009D . Here we provide 
a brief summary. 

Photometric redshifts were computed for every ob- 
ject in the EDisCS fields using two independent codes, 
a mo dified version of the publicly av ailable Hyperz 
code (iBolzonella. Miralles. fc Pelldll2000f) and the code 
of Rudnick ct al. (2001) with the modifications presented 
in lRudnick et al.l (|2003|). The accuracy of both methods 
is a{Sz) « 0.05 - 0.06, where Sz = ""Y+zflT ' 
ting stellar templates to the observed SEDs of stars we 
searched for zeropoint offsets and found no offsets ex- 
cept for a small one in the B-band of CL1353. 0-1137. 
We applied the offset for this one band when perform- 
ing the photometric redshift fits. We established mem- 
bership using a modified versi o n of the technique first 
developed in Brun ner fc Lubii] |2000'), in which P{z) is 
integrated in a slice around the cluster redshift for the 
two codes. The width of the slice around which P{z) is 
integrated should be on the order of the uncertainty in 
redshift for the galaxies in question. In our case we use 
a Az = ±0.1 slice around the spectroscopic redshift of 
the cluster Zdust- We reject a galaxy from our member- 
ship list if Pciust < Pthresh for either code. We calibrate 
Pthresh from our ~ 1900 spectroscopic redshifts. Our val- 
ues of Pthresh Were chosen to maximize the efficiency with 
which we can reject spectroscopic non-members (down 
to / = 22) while retaining at least « 90% of the con- 
firmed cluster members, independent of their rest-frame 
{B — V) color or observed {V — I) color. In practice 
we were able to choose thresholds such that we satisfied 
this criterion while rejecting 45-70% of spectroscopically 
confirmed non-members. Applied to the entire magni- 
tude hmited sample, our thresholds reject 75 — 93% of 
all galaxies with /tot < 24.9. It is worth noting that 
it is very difficult to assess our absolute contamination 
for two reasons. First, even the extensive spectroscopy 
we currently have was performed on a subsample of the 
photometric catalog that was designed to exclude objects 
with an extremely low probability of being at the clus- 
ter redshift. Any estimates based on this spectroscopy 
may therefore not be entirely indicative of the true con- 
tamination down to the spectroscopic completeness limit. 
Second, we do not have spectroscopy for galaxies down to 
the faint limit of the photometric catalog and it becomes 
impossible to definitively measure the contamination at 
these faint magnitudes without significantly deeper spec- 
troscopy or highly model dependent assumptions. 

Our method establishes cluster membership using a 
redshift interval smaller than that employed in other 
photometric-based membership techniques (e.g. Toft et 
al. 2004) and therefore should suffer considerably lower 



field contamination. As a check of how much more con- 
tamination we wo uld have if we adopted the technique of 
iToft et al.l (I200I we have re-measured our membership 
requiring that each galaxy be within Az = ±0.3 of the 
cluster redshift. The number of cluster members with 
this technique is typically 2-3 times larger than when us- 
ing our membership technique, implying a correspond- 
ingly larger contamination. 

Despite the apparently good performance of the pho- 
tometric redshift technique, the Zphot estimates are only 
well tested at relatively bright magnitudes, e.g. / < 22. 
Because the Zp^io^-based membership technique is largely 
untested at / > 22 it will be difficult to trust the faint end 
slope of the LF derived from such techniques. For this 
reason it is desirable to use complementary photometric 
methods to establish membership. 

3.2. Statistical Background Subtraction 

An independent method of establishing cluster mem- 
bership is the statistical subtraction technique (e.g. 
Aragon-Salamanca et al. 1993; Stanford et al. 1998). 
In this technique, number counts in the cluster field are 
compared to those in an "empty" field and the excess 
counts are used to assign a membership probability to 
each galaxy in the cluster field. This method becomes in- 
creasingly inefficient at high redshift, where the contrast 
of the cluster against the background becomes increas- 
ingly low. In addition, this method provides no member- 
ship probability for individual galaxies, but rather gives 
every galaxy in a given region of magnitude (and color) 
space an identical probability. At the same time, it suf- 
fers from completely different uncertainties than the pho- 
tometric redshift technique and is a useful complement 
to judge the robustness of our results. 

Ideally, the comparison catalog used to create the field 
counts should contain the same bands as used in the 
cluster fields and cover a large enough area to mini- 
mize cosmic variance. For our statistical-background- 
subtraction-based membership we utilized a "field" cat- 
alog from the Canada France Deep Field (CFDF; Mc- 
Cracken et al. 2001)^^. This field has the advantages of 
having matched aperture photometry in V and /-bands 
and AUTO magnitudes in the /-filter, while also covering 
0.25 square degrees, roughly 20 times the area of the op- 
tical coverage in an individual EDisCS field. The depth 
of the CFDF is only / — 24.5 and so all LFs computed 
via statistical background subtraction will be limited to 
/ < 24.5. The CFDF is the only publicly available field 
that satisfies our requirements for a background field. 
These were: 1) that it must have photometry in at least 
V and I since these filters are in common for both of 
the EDisCS filter sets (BVIK and VRIJK) and 2) that it 
must have a large enough area to overcome the effects of 
cosmic variance in the background subtraction estimate. 
While there are other fields with deep multi-field pho- 
tometry over a moderate area (e.g. Chandra Deep Field 
South, NO AO Deep Wide Field Survey), there are no 
publicly available surveys with both deep V and I at a 
depth comparable to EDisCS and with large enough area 
to overcome cosmic variance. For example. The Chan- 
dra Deep Field South (CDF-S) that was targeted by the 
FIREWORKS survey (Wuyts et al. 2008) is known to be 

This catalog has been kindly provided to use by H. McCracken 



6 



Red Sequence Cluster Galaxy LF at z < 0.8 



underdense at z ^ 0.7 compared to the much larger Ex- 
tended Chandra Deep Field South (ECDF-S; Taylor et 
al. 2009) and so is not a good sample of the mean back- 
ground. Also, the NOAO deep wide field survey (Brown 
et al. 2007) , which we use in §6.21 has no V filter and a 
very wide _B-band filter (essentially U + B), which makes 
it impossible to use as a background field for the EDisCS 
clusters with only BVIK photometry. 

We use a ni e thod similar to the one presented by 
iPimbblct et al.l 1)200 2') and refer to that paper for de- 
tails, although we summarize it briefly here. We bin the 
CFDF data and our own in observed {V — T) color and 
Iauto using bins of 0.5 in color and magnitude (using 
color bins of 0.3 results in nearly identical LFs). Note 
that we do not use /tot when performing the statisti- 
cal subtraction, as the CFDF does not have aperture 
corrected magnitudes. We assume that the AUTO mag- 
nitudes perform similarly for both surveys. In a given 
bin we scale the number of field galaxies to the area of 
the cluster under consideration to derive the number of 
expected field galaxies. We first retain all spectroscop- 
ically confirmed members and exclude all spectroscopi- 
cally confirmed non-members. Then we subtract off a 
random subset of the remaining galaxies equal in num- 
ber to the expected number of field galaxies (minus the 
number of spectroscopically confirmed non-members) to 
obtain a realization of the cluster member population. 
In bins where the number of expected field galaxies are 
greater than the number of member candidates, we merge 
adjacent bins in color until the number of expected field 
galaxies is greater than or equal to the number of mem- 
ber candidates in the expanded bin. This is analogous to 
expanding the bins until the membership probabilities 
agai n lie between a,nd 1 . As explained in Appendix 
A of iPimbblet et all (|2002f ) this method has an advan- 
tage over similar methods in that it preserves the origi- 
nal probability distribution, albeit smoothed over larger 
scales. 

The moderately large area of the CFDF gives an accu- 
rate representation of the mean density of field galaxies 
but on spatial scales similar to that of our clusters the 
number counts of field galaxies may vary and the true 
underlying field may be systematically different from the 
mean. We use the entire CFDF area to calculate our 
best estimate of the membership sample for each cluster. 
When calculating the uncertainty in the cluster member- 
ship, we split the CFDF into tiles, with each tile having 
the same area as the area of the cluster under considera- 
tion. In practice this resulted in greater than 20 indepen- 
dent tiles in the CFDF. We then performed 100 Monte 
Carlo iterations of the subtraction, where each iteration 
uses a randomly chosen tile to derive the expected field 
population. 

4. MEASURING THE LUMINOSITY FUNCTION 

In this section we will present our niethod for deter- 
mining rest-frame luminosities, for measuring the LF of 
cluster galaxies as a whole and split by color, and for 
fitting Schechter (1976) functions to the measured LFs. 
We will present a comparison of LF^p and LFgs and dis- 
cuss why robust LF determination of cluster galaxies can 
only be made for the red galaxy population. 

4.1. Determining Rest-frame Luminosities 



Rest-frame luminosities L^^'"'* and rest-frame col- 
ors were calcu l ated using the technique described in 
iRudnick et all (|2003D and assuming that every galaxy 
selected as a cluster member has z — Zdust ■ Our L^f ^' es- 
timates were computed from the matched aperture mag- 
nitudes (see §2.2p . which almost certainly miss flux com- 
pared to the /tot estimate. To scale our U^^^ estimates 
to total values we therefore multiply every U^^^ value 
by the ratio of the total /-band flux to that in the /- 
band matched aperture. The median correction ranges 
from a few percent at /tot ^ 20 — 21 to ^ 30 — 50% at 
/tot - 24.4 - 24.9. 

Which rest-frame luminosities we are able to use de- 
pends on which technique we employ to determine clus- 
ter membership. For the photometric redshift method 
the full range of rest-frame wavelengths are available, as 
the probability of each galaxy residing at the cluster red- 
shift is computed directly from its SED. Therefore the 
SED is by definition consistent with being at (or near) 
the cluster redshift and any interpolation between the 
observed bands based on the templates at that redshift 
should yield a robust estimate of U^^*^. We therefore 
can compute rest-frame magnitudes of cluster members 
in many rest-frame bands spanned by our observed filter 
sets, e.g. gresti frest, and irest- The rest- frame NIR lumi- 
nosity functions will be presented in Aragon-Salamanca 
et al. (in preparation). 

The statistical background subtraction method, how- 
ever, limits the rest-frame wavelengths for which lumi- 
nosities can be robustly computed to those that are 
straddled by the observed subtraction filters. The rea- 
soning is as follows. Recall that the photometric red- 
shift technique uses the full SED information to deter- 
mine membership on an individual basis. With statis- 
tical subtraction, however, the membership probability 
is not known for each galaxy, but rather for all galaxies 
in a region of color-magnitude space based on their rel- 
ative numbers with respect to those in an empty field 
image. This implies that some fraction of the galax- 
ies classified as members will actually be at different 
redshifts than the cluster. For rest-frame wavelengths 
straddled by the observed subtraction bands (in our case 
Ay < (1 -|- Zciust) X Arest < A/) this is not a problem, 
as the color of every candidate member is constrained to 
be similar to that of the very cluster galaxies that cause 
the overdensity in counts in that color-magnitude bin, re- 
gardless of whether or not that candidate truly is a mem- 
ber. Therefore the use of templates at Zdust can be used 
to determine L^f^' without large systematic errors if the 
galaxy is truly a non-member. However, this statistical 
subtraction method docs not insist that the SED of the 
galaxy outside of the observed subtraction bands is con- 
sistent with one at the cluster redshift. For this reason, 
rest-frame wavelengths outside of the subtraction bands 
will be subject to uncertain extrapolations and will not 
be robust. For clusters at our redshift, the condition 
Ay < (1 -|- Zdust) X Aresi < A/ is approximately satisfied 
for the Qrest and B^est-bands, which we limit ourselves 
to for LFs computed with statistical subtraction. 

4.2. A Non-parametric Estimate of the LF 

We first measure the LF of every cluster by simply bin- 
ning the sample into absolute magnitudes and counting 
the number of galaxies in each bin. As is done in previous 
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works, we exclude the Brightest Cluster Galaxy (BCG) 
and galaxies brighter than the BCG from the LF compu- 
tation. The properties of the EDisCS BCGs haven been 
presented separately in lWhilev et al.l (|2008f) . 

For LFzp the error bars in each bin represent the Pois- 
son errors o n the re t ained galaxies, computed using the 
formulae of lGehrelsl (|1986D . For LF^s the best-fit LF is 
that derived using the subtraction over the whole CFDF. 
There are two sources of error that contribute to LFgs. 
The first source is the Poisson error on the number of 
galaxies in each cluster field retained as members. The 
second source of error originates in the uncertain back- 
ground measurement, which we determine using Monte 
Carlo realizations for small sub-tiles of the CFDF in esti- 
mating the field (see . In this case we computed the 
LF for each Monte Carlo realization of the subtraction 
and took the 68% confidence intervals of the resultant 
LFs as an estimate of the error. This error was then 
added in quadrature to the Poisson error to achieve a 
total error. 

In constructing the LF for each cluster there are two 
issues to consider, the detection limit in observed total 
magnitudes and the corresponding limit in absolute mag- 
nitude. As described in W05 we establish our complete- 
ness in observed Lband magnitude in an empirical way by 
comparing our number counts to those from much deeper 
surveys (see W05, Fig. 1). There is ample evidence that 
the intrinsic slope of the I-band number counts is a ris- 
ing power law at faint magnitudes (e.g. Metcalfe et al. 
2001; Heidt et al. 2003) and we define our completeness 
as the magnitude at which our observed number counts 
in total magnitudes deviate from a power law defined by 
the deeper observations. There are two reasons this is 
reasonable. First, the number counts contributed by the 
cluster at faint magnitudes is much smaller than the con- 
tribution by the field. This is evidenced by the fact that 
80—90% of galaxi es are rejected b y statistical subtraction 
at /tot < 24.9 (jPello et al.ll2009f) . Also, the slope of our 
number counts is parallel to that from deeper fields at 
22 < /tot < 24 for the high-z clusters and 23 < /tot < 24 
for the low-z clusters, where we expect the cluster to no 
longer contribute significantly to the counts. For this 
reason we feel that our faint counts can be directly com- 
pared to that of the field. Second, our total magnitudes 
(which include an aperture correction) result in a rapid 
drop off in the number counts at faint magnitudes. This 
is not seen in surveys that measure magnitudes with- 
out an aperture correction but is a direct result that 
we count for a mini mal amount o f mis sing flux in our 
faintest galaxies (see 'Labbe et al.' (2003|) for a more de- 
tailed explanation) . Labbc et al. ( 2003) also showed that 
a limit defined in this way corresponds to a near perfect 
detection probability. Because this is a rather conser- 
vative estimate of our completeness the signal-to-noise 
is still high (typically > 10; W05) all the way down to 
our detection limit, allowing the robust computation of 
magnitudes and colors. 

Once we have established our completeness limit in ob- 
served magnitude we translate this, for every rest-frame 
filter, into an absolute magnitude limit that is the most 
conservative (i.e. brightest) given the whole range of pos- 
sible galaxy SEDs. If a redshifted rest-frame filter for a 
given cluster redshift is blueward of the observed /-band 
the brightest limit corresponds to that computed using 



a 10 Myr old si n gle ag e population with solar metallicity 
and a lSalpete"il ()1955h IMF. This is perhaps overly con- 
servative for red galaxies, but results in the most conser- 
vative limit for the whole catalog, so that we are equally 
complete at all galaxy colors. For redshifted rest-frame 
filters redward of the observed /-band we used E lliptical 
template from lColeman. Wu. fc WeedmanI ()1980[ ). 

We also created composite LFs for subsamples split 
by redshift and cluster velocity dispersion. We cre- 
ated th e co mposite and its error using the method of 
CollessI (figsai . which was also discussed in detail in 

mi 



Popesso et al.l ()2005[ ). With this method, the composite 



LF at every magnitude represents the mean fraction of 
galaxies compared to the number in a normalization re- 
gion. We choose the normalization region to be all mag- 
nitudes brighter than the brightest completeness limit 
that all clusters in that subsample have in common. 

When creating the composite clusters, we correct them 
for passive evolution to the mean redshift for that sub- 
sample. As we will describe in subsequent sections, only 
the LF for red cluster galaxies can be robustly deter- 
mined and we concentrate mostly on those for the rest 
of the paper. DL07 showed that the colors of the red 
sequence can be well fit by a passively evolving model 
with Zform '^2 — 3. We correct the rest-frame magni- 
tudes using a Zform — 2 single stellar population (SSP) 
Bruzual & Chariot (2003; hereafter BC03) model with 
Z = 2.5Zq. In practice, this small evolution correction 
does not change the binned LF with respect to that com- 
puted with no correction. This is because the amount of 
evolution from each cluster to the center of its redshift 
bin is significantly smaller than the 0.5 magnitude bin 
size used in constructing the LF. For the same reason 
the exact choice of model used makes little difference in 
the resulting composite LF. 

We compute the LF in two different physical radii, r < 
0.75Mpc and r < 0.5/?2oo, where /?200 is defined as the 
radius within which the density is 200 times the critical 
density: 



R200 = 1.73 



1000 km s-i ^orrTvmF 



^100 Mpc 



wh ere a is the cluster velocity dispersion (|Finn et al.l 
[2005h . The area defined by these two radii is entirely 
contained within the EDisCS fields for all but one of 
our clusters (CL1227.9-1138) for which we take only the 
inscribed area into account when performing the statis- 
tical subtraction^". For this cluster the lack of data for 
~ 50% of the galaxies within 0.5/?2oo and ~ 60% within 
0.75Mpc should not bias the values of M* but will re- 
sult in a larger error bar on that value. For only two of 
the most massive clusters, CL1216. 8-1201 and CL1232.5- 
1250, is 0.5/?2oo larger than 0.75Mpc. Our conclusions 
are insensitive to the exact choice of radii and unless 
otherwise stated we will use r < 0.75Mpc as it is most 
always the larger of the two and hence will produce the 
highest S/N LF. 

20 Using the EDisCS data it was realized that the LCDCS BCG 
candidate for CL1227.9-1138 was not the actual BCG. The true 
BCG is located significantly ofTcenter in our FORS data, resulting 
in the loss of area. 
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4.3. Schechter Function Fits 

We fit iSchechteil (119761 ) functions to tlic binned LFs in 
each cluster. To fit we created a coarse grid in the three 
fitted parameters, i.e. cj)* , M* , and a. We calculated the 
value at each grid point and took the best-fit solution 
as an initial guess for the parameters. We then re-fit the 
parameters with a narrower range and a finer sampling 
in the parameter space. We determined the formal un- 
certainty on each parameter by first converting the at 

2 Icy 

each grid point into a probability via P,p*.M*.a — e^-^ ' 
and then by marginalizing the probability along the other 
two parameters to obtain a probability distribution for 
the parameter in question. We then measured the limits 
in this parameter that enclosed 68% of the probability as 
the 1-sigma formal error bar. 

To assess the reliability of such fits we created a set of 
mock binned luminosity functions by randomly drawing 
from a set of input values, i.e. the number of galaxies, 
M*, and a. The errors on each mock LF were Poisson er- 
rors on the number of galaxies in each bin. For a given set 
of parameters we created 100 mock realizations of that 
LF and fit each realization using the procedure above, 
and over the absolute magnitude range present in our 
data. While all three Schechter parameters are highly 
degenerate, we found that the most poorly constrained 
parameter was a followed by M* . The ability to retrieve 
the parameters was also dependent on the input value 
of a, since steeper (more negative) a's produced more 
biased answers. For the red galaxies to which we limit 
our analyses (see subsequent sections) a > —0.6 and the 
bias produced by a steep slope is not severe. Nonethe- 
less, through these simulations we found that it was im- 
possible to constrain all three parameters simultaneously 
using the data from an individual cluster, or even from 
a composite luminosity function of only a few clusters. 
We did find however, that we could constrain all three 
simultaneously if we fit a LF with characteristics akin 
to the composite LF of the entire EDisCS sample, split 
into two bins of redshift. We therefore derive a and its 
uncertainty for the entire EDisCS sample for each band 
in each redshift bin and use that a when fitting the in- 
dividual and stacked luminosity functions when split by 
velocity dispersion. Even when fitting to the whole sam- 
ple, however, the uncertainties on a are non-negligible. 
To account for this uncertainty in the fitting of individual 
clusters or subsamples of the EDisCS clusters, we fit the 
Schechter function to the data 100 times, with a fixed 
each time but drawn randomly from a Gaussian with a 
mean and sigma taken from the fit to the total stacked 
cluster sample. The 68% confidence interval in the dis- 
tribution of M* from these 100 iterations was then added 
in quadrature to the formal uncertainties, derived with a 
fixed a, to derive the total uncertainty in M* . This may 
overestimate the error in M* as it includes some of the 
sampling error twice. 

4.4. Splitting LFs by color 

We divide our sample by {V — /) color into red se- 
quence galaxies and bluer galaxies. For each cluster we 
fit the zeropoint of the color-magnitude relation (CMR) 
in {V — I) assuming a fixed slope of - 0.09 and using the 
outlier resistant Biweight estimator (jBeers "eFall [l990l) 
for the zeropoint. In performing the fit we only use spec- 



troscopically confirmed cluster members with no emis- 
sion lines. This was the same method as used by DL07. 
We give the best-fit zeropoints in Table [T] for the 16 clus- 
ters for which a robust LF determination is possible (see 
iJS])^^. A relatively constant slope of the CMR can be un- 
derstood if the slope is primarily a result of a metallicity 
trend with magnitude (e.g. Kodama & A rimoto 1997) 
among galaxies with a uniformly old age (|Bower et al.l 
1992), at least among br ight gal axies. As sh o wn in , e.g. 
iKodama fc Arimotol ([19971 and Bower et al.l (|1998f ). the 
rate of change of color with time is insensitive to metal- 
licity, so using the local value for the CMR slope with 
our intermediate redshift clusters is a reasonable assump- 
tion. As in DL07, we select red sequence galaxies as those 
within ±0.3 magnitudes of the best-fit CMR. This is a 
compromise between the completeness and purity of our 
red sequence sample. By allowing our color cut to ex- 
tend below the CMR, we ensure that we do not miss 
red galaxies that are slightly bluer than the CMR, but 
also increase the possibility that there may be some blue 
galaxy contamination at fainter magnitudes where our 
photometric errors increase. We tested the sensitivity of 
our results were to the exact form of our red sequence 
select in two ways. First, we varied the width of our se- 
lection slice by ±0.05 mag. This corresponds to the « 0.1 
magnitude error in (V — I) for ga laxies at the EDisCS 
magnitude limit (jWhite et al.ll2005h . Second, we selected 
all galaxies redward of the CMR and then reflected them 
across the CMR. This latter m ethod is similar to what 
is used bv lGilbank et al.l (|2008f ) and insures high sample 
purity at the risk of missing intrinsically bluer/younger 
galaxies still formally on the red sequence. In all cases 
the we find that the LFs with these different methods 
are consistent to within 1-sigma, indicating that our re- 
sults are robust against variations in the red sequence 
selection. We believe that this must be partly true due 
to our conservative magnitude limit and extremely deep 
VLT photometry. 

For each of the samples split by color we compute the 
individual and composite LFs as described above. As 
shown in DL04 and DL07, it is also important when es- 
tablishing the effective magnitude limit on the red se- 
quence to take into account that the S/N of the color 
measurement of galaxies becomes worse for redder galax- 
ies at a fixed /tot (see Fig. 1 of DL07). We take this into 
account when determining our completeness limit and 
find that we may be missing some red sequence galaxies 
in the 24.4 < I < 24.9 magnitude bin. Although our 
LFzp estimates for the high-z clusters are computed to 
/ = 24.9, all of the trends described in this paper are 
completely dominated by effects in the bins at / < 24.4. 
We therefore do not worry about this minor incomplete- 
ness in our last bin. 

4.5. A Comparison Between Methods 

We assess the robustness of our LFs by comparing LF^p 
and LEgg. In Figure [5] we show the composite LF of 
all EDisCS clusters as computed with the two methods. 
The LFss of all galaxies has a steeper faint end slope and 
an overabundance of bright galaxies compared to LF^p . 

Our values are given at /tot= whereas those from DL07 were 
given at an apparent magnitude that corresponds to My= —20 
when evolution corrected to z = 0. DL07 also use Vega magnitudes. 
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This same behavior is apparent, albeit at lower signifi- 
cance, in all the composite and individual LFs. We also 
compute the LFs separately for blue and red galaxies and 
plot these in the middle and right panels of Fig [2] respec- 
tively. It is obvious from this figure that the discrepancy 
only exists for the blue galaxies. In contrast, LFgs and 
LFzp agree completely for red galaxies, as was found in 
DL07. 

There are at least two possible reasons for the large dif- 
ference in the faint end slope between the two techniques 
that only manifests itself for blue galaxies. First, the ef- 
fectiveness of LFss is critically dependent on the validity 
of the field counts used to make the statistical subtrac- 
tion. The faint-end slope of the blue number counts is in 
general steep (e.g. Koo 1986) and we have checked that 
the faint-end slope of the counts in the CFDF is signifi- 
cantly steeper for blue than for red galaxies. Because the 
faint-end slope of the blue galaxy counts is so steep, the 
faint-end slope of the cluster LF is critically dependent 
on the exact value of the slope. Specifically, the faint- 
end slope of the counts in the comparison field needs to 
be the same as the faint-end slope of the counts for field 
galaxies in the cluster field. If there are slight differences 
in the way that magnitudes are measured between the 
field and cluster catalogs, an incorrect faint-end cluster 
LF can be measured. Indeed, although AUTO magni- 
tudes are used for both the cluster and CFDF catalogs 
the seeing of the CFDF catalog is ~ 1.5 — 2 times worse 
than that of the EDisCS catalogs and there has been 
no attempt to match SExtractor catalog parameters. As 
a result, magnitude dependent differences in the AUTO 
magnitudes could be present between the two catalogs 
and this could cause the very steep faint end slope of 
LFss for blue galaxies. We have checked that a magni- 
tude independent change in the CFDF magnitudes of up 
to 0.2 magnitudes has no appreciable effect on the faint 
end slope but have not explored more complicated mag- 
nitude dependent effects. We conclude that differences 
in the way the two surveys measure magnitudes makes it 
difficult to measure the faint-end LFgs for blue galaxies, 
where the magnitude measurement of faint galaxies is so 
crucial. The faint-end slope of red field galaxies in the 
CFDF is much shallower and small errors in the magni- 
tude measurements will not lead to as large of errors in 
the LF. 

Another possible reason for the difference between LFgs 
and LFzp for blue galaxies, specifically the large differ- 
ence in the faint end slope, may come from limitations 
in the photometric redshift techniques. For the spec- 
troscopic sample we verified that the photometric red- 
shifts performed similarly for red and blue galaxies. Un- 
fortunately given the spectroscopic magnitude limit we 
were not able to verify how the photometric redshifts 
performed at faint magnitudes. In general, the perfor- 
mance of photometric redshift codes depends on the S /N 
of the flux measurements since a higher S/N measure- 
ment allows for a better localization of the features (e.g. 
the 4000 A break) used to determine the redshift. For 
galaxies with weaker intrinsic features in their SEDs, e.g. 
blue galaxies, the photometric S/N must be higher to 
yield a comparable redshift accuracy as for galaxies with 
stronger features, e.g. red galaxies with strong 40G0A 
breaks. Since we determine cluster membership by in- 



tegrating P{z), a poorer constraint on Zphot with a cor- 
respondingly broader P{z) will result in a Pdust that 
may fall below the Pthresh value that was calibrated for 
brighter galaxies. As an additional complication, the 
slope of the blue star-forming galaxy sequence (the blue 
"cloud") is such that faint blue galaxies are typically 
bluer than bright blue galaxies, meaning that the photo- 
metric redshifts will perform correspondingly worse. To 
assess whether this effect could cause the downturn on 
the faint-end LF^p for blue galaxies we examined the de- 
pendence of the Zphot 68% confidence limits on Mg for 
blue and red galaxies with Zphot = Zdust ± 0.05. For 
red galaxies the internal Zphot errors remain small and 
increase only slowly with Mg. For blue galaxies, how- 
ever, the internal errors rise more rapidly with increas- 
ing Mg and there is a population of blue galaxies with 
very large errors. Both the blue galaxies with very large 
errors and those on the upper envelope of the main error- 
magnitude relation are flagged as interlopers by the pho- 
tometric redshifts. The absolute magnitude where this 
increase in the Zphot uncertainties of blue galaxies oc- 
curs coincides with the magnitude where the faint end 
slopes of LFss and LF^p start to diverge. The difficulty 
in using Zphot to establish membership at fain t magni- 
tudes is explored further in iPello et al.l (|2009| ). It may 
be that the best way to study blue galaxies with photo- 
metric techniques is by using a combination of statistical 
background subtraction and photometric redshift mem- 
bership techniques, such that the photometric redshifts 
are used as a first-pass membership method and the sta- 
tistical background subtraction is then used to subtract 
off any residual (e.g. Kodama et al. 2001; Tanaka et al. 
2005). In practice this will require either a large field 
sample with identical photometry (and hence photomet- 
ric redshift performance) as the target field or a cluster 
image with a wide enough area to have minimal contam- 
ination from the cluster at the outskirts of the image. 

As mentioned, these two problems in determining the 
faint end should not be (and apparently are not) as se- 
vere for red galaxies as for blue. Photometric redshifts 
seem to perform better for red galaxies than for blue, at 
least in the realm of decreasing photometric S/N. The 
source of the discrepancy between LFgs and LF^p at the 
bright end is not as clear. The CFDF appears to be 
slightly underdens e with respect to the FORS Deep Field 
(jHeidt et al.||2003|) and the COMBO- 17 number count s 
from The Chandra Deep Field South (|Woff et al.l l2004V 
which would serve to increase the LFgs value for EDisCS. 
Also, despite our best efforts at calibration of Zphot for 
bright sources from the spectroscopic sample, the pho- 
tometric redshifts may reject a slightly larger number of 
blue members than red members, which would push LFzp 
down. In the end we must conclude that the determina- 
tion of the blue-galaxy cluster LF is not robust when only 
using photometric redshifts or statistical subtraction. 

The red galaxy LFs, however, agree astonishingly well, 
indicating that the red galaxy LF is robust to the exact 
method used. We therefore limit most of our subsequent 
analysis to the red galaxies only. 

4.6. The Local Luminosity Function 

To measure evolution in the LF it is important to 
have an appropriate local sample. For many parame- 
ters of the galaxy population, e.g. the star-forming frac- 
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Fig. 2. — The grest-band composite luminosity functions of EDisCS cluster galaxies. The left panel is for all galaxies, regardless of color. 
The middle panel is for blue galaxies and the right panel is for red galaxies (see text for definition of colors). The squares show the LF 
determined using statistical subtraction LFss and the triangles show the LF determined using photometric redshifts LFzp. The solid and 
dotted curves show the best-fit Schechter function fits to LFzp and LFss respectively and the vertical arrows of the same line type show 
the corresponding best-fit values of M* . The horizontal error bars at the base of the arrows give the 68% confidence limits in M* . When 
including all galaxies LFss has a steeper faint end slope and a larger number of bright galaxies than LFzp • These difference can be traced 
to the blue galaxies. Both techniques give identical results for the red galaxies. 



tion (Poggianti et al.l l2006[ ) and the early-type fraction 
([Pesai ct al. 2007), there is a strong dependence on a at 
intermediate and high redshift, implying that evolution 
can only be measured in samples matched in velocity dis- 
persion. No dependence of the LF of al l cluster galaxies 
on a has been found at low redshift (|de Propris et al.l 
|2003| ) but we wish to test this for red sequence galax- 
ies specifically at intermediate and high redshift. For 
our purposes we therefore require a local sample that 
has the same range in cr as our sample and allows for 
the computation of a luminosity function just for red se- 
quence galaxies. It is also desirable that enough local 
clusters be used so as to average over cluster-to-cluster 
variations and minimize the uncertainties in the local 
anchor of any evolutionary trends. Finally, it is advan- 
tageous if the local LF has been computed in multiple 
bands, to allow the measurement of wa v elength depen- 
dent e voluti on. Ide Propris et al.l (|2003l ). iPopesso et al.l 
(|2005D . and iPopesso et alT i 20061 ) computed composite, 
high S/N LFs from the 2dFGRS and SDSS respectively. 
Ide Propris et al.l (|2003f l compute their LFs only in the bj- 
band and do not compute them as a function of galaxy 
color. IPopesso et al.l ()2006[) presented composite LFs for 
X-ray selected clusters in multiple bands and as a func- 
tion of galaxy color; however, we choose to construct our 
own SDSS luminosity function, for the following reasons. 
The sample of (Popesso et al. 2005;2006) is x-ray se- 
lected, which may cause biases in the comparison of the 



local sample to t he EDisCS sarnple, w hich is optically 
selected. Second, IPopesso et al.l ()2006D split their LFs 
by color, but not in an analogous way to the EDisCS 
sample, which again complicates t he comparison to our 
results. Finally, the raw LFs from IPopesso et ahl ()2006D 
are not published, but only the two-component Schechter 
fits, which also complicates the comparison to our LFs. 

Our cluster sample is a subs et of the sample presented 
in Ivon der Linden et al.l (|2007l ) . This parent sample was 
selected from the C4 catalog of lMiller et"all (j2005), but 
employs improved algorithms to identify the BCG and 
measure the velocity dispersion. We limit our analysis 
to clusters at 2; < 0.06, to ensure that the individual 
cluster LFs are complete down to the passively evolved 
hmit of the EDisCS clusters (see below), which results 
in a sample of 167 clusters. With this redshift cut-off we 
can limit our analysis to galaxies with r < 20, where the 
star/galaxy separation is still robust and where colors 
can be robustly determined. We used a global field sam- 
ple drawn from the SDSS DR4 catalog and use the model 
magnitudes to measure colors and cmodel magnitudes 
to measure the total magnitude. Colors are measured 
with the model magnitude, since that measure adopts 
the same aperture in different bands. The cmodel mag- 
nitudes fit a linear combination of an exponential and 
de Vaucouleurs profile to each galaxy and integrate the 
combination of these two to derive a total magnitude. 
It is well known that the Petrosian magnitude of SDSS 
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misses flux, especially for early type galaxies with de 
Vaucouleurs profiles and cmodel magnitudes should be 
closer to total. 

We isolated cluster members using a statistical sub- 
traction technique similar in principle to what was used 
for the EDisCS clusters, but with some significant modi- 
fications. Since photometry of identical depth and band- 
pass coverage was available for both our local cluster 
and "field" samples in multiple bands, we performed our 
statistical subtraction in 4 dimensions, using bins in r- 
magnitude as well as bins in g — r, r — i, and i — z, and 
using 0.5 bins in magnitude and 0.3 magnitude bins in 
color. This technique has two main advantages over the 
2D, i.e. magnitude and single color, subtraction used in 
EDisCS. Although interloper galaxies may have identi- 
cal colors to cluster members in one color, as more colors 
are considered, it is increasingly difficult to mimic the 
colors of galaxies at the cluster redshift, thereby increas- 
ing the contrast of the cluster against the background 
and reducing the number of contaminating galaxies. In 
addition, because we require that retained galaxies have 
colors matching that of cluster members from g all the 
way to z, it is possible to robustly determine rest-frame 
magnitudes for any bands in betwe e n, e.g . grest, rrest, 
and irest- Following 'Pimbb let et al.l (|2002D here we also 
expand the bins in magnitude first and then in color if 
the number of expected field galaxies exceeds the number 
of candidate members. We determine rest -frame magni- 
tudes using the kcorrect software v4.1.4 (|Blanton et all 
|2003| ). For each cluster we performed 100 Monte Carlo 
realizations of the subtraction, where each iteration used 
a new set of random numbers which were then compared 
against the membership probabilities to determine clus- 
ter membership. The distribution of the LFs for the full 
Monte Carlo simulation was used measure the uncertain- 
ties in the LF. 

To isolate the red sequence in the SDSS clusters we per- 
formed an outlier-resistant fit the to CMR of the compos- 
ite SDSS cluster population using rest-frame magnitudes 
and colors and only fitting spectroscopically confirmed 
members with no Ha emission^^. Evolution corrections 
to the mean redshift for red galaxies are very small over 
our redshift range and do not make a difference when 
fitting the CMR or deriving the LF. Due to the large 
number of galaxies over a large range in absolute mag- 
nitude, we were able to simultaneously fit the slope and 
zeropoint of the relation. For measuring the greet and 
rrest LFs we fit the CMR in (g — r)rest vs. grest and 
{g — r)rest vs. rrest respectively. For measuring the irest 
LF we fit the CMR in (r — i)rest vs. irest- Similar as to 
what was done for the EDisCS sample, we then classi- 
fied as red galaxies every galaxy within a stripe centered 
on the CMR in the color and magnitude used above for 
each rest-frame band. The width of this stripe in each 
color was chosen to correspond to A{U — V)rest = 0.3 
at z = 0.6 assuming that the scatter in the CMR is due 
entirely to age. 

We created individual and composite SDSS cluster LFs 
using identical procedures as with the EDisCS clusters. 

The presence of Ha was indicated by a measurement of 
Ha with S/N(Ha)> 5. The Ha measurements were obtained 
from http : //ww w . mpa-KarchinR ■ mpg ■ de/S DSS/DR4/ raw.data ■ html 
as computed by IBrinchmann et al.l I I2004I) 



i.e. binning the individual clusters in a bsolute magni- 
tude and creating the composite following lCollessI (|1989f) . 
For every Monte Carlo realization of the subtraction we 
computed the Poisson uncertainty on the LF. As with 
EDisCS the LFs were computed within r < 0.75Mpc 
and r < 0.5i?2oo (as computed using Eq. [T]). To deter- 
mine the absolute magnitude limit down to which we 
construct the LF while probing the same galaxies as in 
the EDisCS clusters, we took the absolute magnitude 
limit for our highest redshift EDisCS cluster and cor- 
rected it for passive evolution down to z = 0.06 using a 
Z ^ 2.5Zq BC03 SSP model with Zform = 2. For exam- 
ple, the g^est-band limit of -18.5 for the highest redshift 
EDisCS cluster corresponds to a ^rest-band limit of -17.5 
for the SDSS clusters. The r < 20 apparent magnitude 
selection is deep enough such that the absolute magni- 
tude limit is the more restrictive cut for all of our SDSS 
clusters. 

We also split our SDSS clusters by a to match the 
velocity dispersion bins in the EDisCS sample. We de- 
fined the SDSS a threshold taking into account the red- 
shift evolution in a that is predicted from the growth 
of the cluster dark matter halos over tim e. Using the 
expected mass accr e tion history of halos (iBoweil 119911 : 
|LaceiL& Cole"1993),'Poggianti_eLaIl('2006') used the re- 
sults of Wechsler ct al. (200j) and [Bullock et al. (200ll) 
to show that clusters with a = 600 km/s at 2 = 0.6 
will grow into clusters with a = 700 km /s hy z — (see 
Fig. 19 of Mi lvang- Jensen et al.l ()2008f )). We therefore 
divide our SDSS sample into high and low velocity dis- 
persion bins using a — 700kms^^ as the divider. There 
are 159 clusters in the low a bin and 8 in the high a bin. 
The mean and median uncertainty in the SDSS velocity 
dispersions is ±62 and ±55 km/s respectively. 

We fit each Monte Carlo realization of the composite 
SDSS LF with a Schechter function, allowing all three 
parameters to vary. The best-fit parameters come from 
the mean of the realizations and the uncertainties on 
the parameters come from the 68% confidence intervals 
of distribution from all of the Monte Carlo realizations 
summed in quadrature with the mean formal fit errors. 
In Figure[3]we show the SDSS LFs for red sequence galax- 
ies and the corresponding Schechter function fits in each 
band for the whole cluster sample. The Schechter fit 
parameters are also given in Table [21 

5. RESULTS 

In this section we present the LFs for red-sequence 
galaxies in the EDisCS clusters. As we discussed tj4.5l we 
will limit our analysis to the red sequence LF. Through- 
out we will use LF^p since it gives us access to a larger 
range of rest-frame wavelengths and allows us to go 
0.5 magnitudes deeper in our high redshift clusters (see 
HA.l^ . We will first present the individual cluster LFs 
and then the composite LFs split by redshift and veloc- 
ity dispersion. In all cases we do not show results for 
CL1119. 3-1129 since this cluster has no NIR data and 
hence has poorly constrained photometric redshifts. We 
exclude the CL1103.7-1245a and CL1103.7-1245b clus- 
ters aX z — 0.70 and 0.63 since these clusters overlap 
on the sky and are too close in redshift to be decom- 
posed with the photometric redshifts. We also exclude 
CL1103. 7-1245 since its redshift (z = 0.96) is too high 
for our imaging to probe far enough down the LF. 
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Fig. 3. — The rest-frame optical composite luminosity func- 
tions of red sequence galaxies for 167 clusters at z < 0.06 from 
the SDSS. The data points, error bars, and solid line represent one 
of the Monte Carlo iterations of the background subtraction, the 
Poisson errors, and the corresponding Schechter function fit. The 
parameters listed are the mean over the full set of Monte Carlo 
realizations and the error bars on the parameters are the quadra- 
ture sum of the Poisson errors and those from our Monte Carlo 
simulation of the background subtraction errors. 

5.1. Individual Cluster LFs 

In Figure|3]we present the red-sequence LF^p for the re- 
maining 16 EDisCS clusters. The LF of these clusters in 
the (7, r, and i-bands is given in the appendix. In all cases 
we fix the faint-end slope to the value determined from 
the composite LF in the relevant redshift range. In gen- 
eral the Schechter function fits are good with only a few 
clusters having measured LFs that are of too poor quality 
to obtain a reasonable fit, e.g. CL1227. 9-1138. In most 
of the remaining cases it appears that the slope deter- 
mined from the composite LF is an acceptable fit to the 
individual clusters, indicating that a universal LF for red- 
sequence galaxies is possibly in place at these redshifts. 
There are, however, some exceptions, e.g. CL1301. 7-1139 
and CL1037. 9-1243, where the composite faint-end slope 
does not seem to adequately represent the cluster LF. 
There is no significant observed trend of M* with cluster 
redshift. This is in contradiction to the simple expec- 
tation that red sequence galaxies will be brighter in the 
past due to passive evolution. This will be addressed in 
the coming sections. 

5.2. Composite Cluster LFs 

As shown in FigurelH the signal-to-noise of the individ- 
ual LFs arc too low to make any conclusions about trends 
with redshift or velocity dispersion. For this reason we 
create composite clusters splitting the sample into two 
bins ai z — 0.6 and in each redshift bin into two bins of 
velocity dispersion at cr = 600kms^^. In all cases we have 
used the faint-end slope as determined from the compos- 
ite in the same redshift bin. The composite LFs are given 
in Table [3] and the Schechter function parameters for all 



of the composite LFs are give in Table [H 

5.2.1. redshift evolution 

In the top panel of Figure [5] we show the composite 
LFs split by redshift at z = 0.6. All LFs are only plotted 
as faint as they are complete and, for plotting purposes, 
are all normalized to have the same total integrated lu- 
minosity. In the bottom panel we have corrected the LFs 
for the passive evolution expected for a population that 
formed all of its stars at z = 2. This formation redshift 
provides a good fit to the color evolution of the bright red 
sequence galaxies in the EDisCS clusters (DL07). At the 
bright end we do not have high enough signal-to-noise 
to constrain the detailed evolution precisely. Within the 
considerable uncertainties, however, it appears that the 
LFs at the bright end are all similar after applying the 
correction for passive evolution. This indicates that the 
bright galaxy population was mostly in place in clusters 
at z < 0.8. 

In contrast to the bright end, there is dramatic evo- 
lution in the faint-end slope of the LF from z — 0.8 to 
z = 0, in the sense that the number of faint galaxies per 
luminous galaxy are increasing with decreasing redshift. 
Other works over recent years have found similar results 
(DL04; Kodama et al. 2004; Goto et al. 2005; Tanaka 
et al. 2005; DL07; Stott et al. 2007; T anaka et al. 2007; 
Gilbank et al. 2008), althoug h lAn^ on (200 d) finds no 
such result. IStott et all (|2007h and [Gilbank et al] (|2008f ) 
have large samples of clusters but we measure the LF 
to one magnitude fainter than they do at all redshifts. 
DL04 and DL07 use the same data as in this analysis, but 
concentrate on the evolution in Mum/A^faint, as defined 
using a passively evolving luminosity threshold. Our re- 
sults using an independent analysis are consistent at the 
1-sigma level with those from DL04 and DL07. 

We also fit Schechter functions to the LFs at different 
redshifts and plot the resulting fits in the top panel of 
Figure m As seen in Table |4] within the EDisCS sample 
there is strong evolution in a at the 0.5 — 0.7 level at 
a significance of 4 — 5 sigma (depending on the filter). 
The strong evolution continues down to the SDSS sample 
with even higher significance. At the same time there is 
no evidence for evolution in M* with redshift. At first 
glance this seems strange since the bright galaxies are 
already in place in the clusters and appear to be brighter 
in the past. However, the evolution in M* only tracks 
the luminosity evolution of the whole galaxy population 
if all the galaxies evolve in the same way. Given the 
rapid evolution at the faint end, this is definitely not the 
case. The highly degenerate nature of a and M* make it 
hard to interpret them simply as evolution of any given 
part of the galaxy population. As a side note, the lack 
of evolution in M* that we observ e is also consistent a t 
the 1 — 2a level with the results of iTanaka et al.l (|2005f ). 
who also find a strongly evolving faint end. 

Stronger constraints on the joint evolution of the 
Schechter function parameters would be enabled most 
effectively by higher signal-to-noise at the bright end, as 
our sample at fainter magnitudes is statistically quite ro- 
bust. The small number of bright galaxies is a direct con- 
sequence of our limited number of clusters and of their 
modest mass (or richness) . The best way to improve the 
constraints on the Schechter function parameters will be 
to construct LFs for much larger samples of clusters. 
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Fig. 4. — The grest-band luminosity functions of red sequence galaxies in the 16 EDisCS clusters for which LFzp could be determined. 
The clusters increase in redshift to the right and down. The solid curve gives the best-fit Schechter function with a slope fixed to the values 
fit to the EDisCS composite LF in the corresponding redshift bin. The arrow and associated error bar indicates the fitted value of Af*and 
its 68% confidence limits. 

In i jG.li we will discuss this further and present implica- 
tions for the formation histories of red sequence galaxies 
in clusters. 



5.2.2. sample split by velocity dispersion 

We also examined the differences between the red se- 
quence LFs when we split our sample by cluster velocity 
dispersion at cr = 600 km/s. The SDSS sample was split 
at a — 700 km/s corresponding to the expected growth of 
clusters over time (see ^^4.61 ) To assess the differences we 
take the ratio of the LFs in different a bins and show the 



ratio in Figure [51 In this figure trends of the ratio with 
magnitude illustrate differences between the high and low 
velocity dispersion clusters. There are no major trends 
with magnitude but in the 0.42 < z < 0.6 redshift bin the 
low velocity dispersion clusters have systematically more 
faint galaxies relative to bright galaxies, than high ve- 
locity dispersion clusters. Still, formally all of these rela- 
tions are statistically consistent with a co nstant. As part 
of the ir analysis of iVium/-^faint DL07 and lGilbank et al] 
()2008f l also divided their samples by velocity dispersion 
and cluster richness respectively. For reference we draw 
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Fig. 5. — The composite rest-frame g, r, and i LFs of red sequence EDisCS and SDSS cluster galaxies split into bins of redshift. The 
composites at 0.6 < 2 < 0.8 and 0.4 < z < 0.6 have 6 and 10 clusters respectively. All plotted LFs and their fits have been scaled so 
that the total luminosity density of the best-fits are equal. Top row: The LFs are presented without any evolution correction. Vertical 
arrows and horizontal error bars give the value of Af* and the 68% confidence interval. The top, middle, and bottom arrow correspond 
to the SDSS, intermediate redshift EDisCS clusters, and high redshift EDisCS clusters respectively. The points for the two EDisCS LFs 
have been offset slightly in magnitude for plotting clarity. Bottom row: The LFs have all been corrected for passive evolution to the mean 
redshift of the SDSS clusters. All symbols are as in the top panel. The luminous red sequence galaxies are all in place out to z = 0.8 but 
the fainter red sequence galaxies built up dramatically over time. 



vertical lines in the bottom panel of Figure [5] that corre- 
spond to the dividing lines between luminous and faint 
galaxies from DL07. Calculating iViumZ-^faint as in DL07 
we find entirely consistent results for the EDisCS sample 
and plot them in Figure [7l For SDSS clusters we find no 
significant difference between the LFs for clusters of dif- 
ferent velocity dispersions, whereas DL07 found a small 
difference^"^. Our SDSS cluster sample is about twice as 
large as that of DL07, uses a more sophisticated method 
for the background subtraction, and uses model magni- 
tudes instead of Petrosian magnitudes. Because of the 
improvements to our SDSS LF with respect that in DL07 
we believe our result although the conclusions are in any 

23 DL07 split their SDSS sample at 600 km/s but found no 
difference when splitting at 700 km/s, although they only had four 
clusters in their a > 700 km/s bin. 



case not sensitive to the difference. In addition the lack 
of a dependence of iVium/. ^faint on cr for local clu sters 
agrees with the results from lde Propris et al.l ([2Q03) who 
found that their bj LFs were consistent at the 2a level 
for clusters with a greater and less than 800 km/s. It 
must be noted, however, that we are only considering the 
red sequence cluster LF and therefore it is not possible 
to dir ectly compare with the results of Ide Propris et all 
()2003l ). Taken together, our results imply that the higher 
velocity dispersion clusters evolved at a faster rate than 
the low velocity dispersion clusters. A depende nce on 
cluste r velocity dispersion is also seen bvlPog gianti et al.l 
()2006f l who find that the fraction of star-forming mem- 
bers in clusters evolves most quickly in clusters with ve- 
locity dispersions similar to those for the most massive 
50% of the EDisCS clusters. 



Rudnick et al. 



15 



iGilbank et al.1 (|2008D analyzed the dependence of 
Mum/iVfaint On cluster richness and also found that richer 
clusters evolve quicker than clusters with low richness. 
They also find that low-richness clusters at 0.4 < z < 0.6 
have a higher Mum/-^faint than high-richness clusters. 
This is in the opposite sense of the trend found by DL07 
and our work, in that higher velocity dispersion clusters 
have a higher iVium/iVfaint than low velocity dispersion 
clusters. While this difference in our work is not very sig- 
nificant, the d i fferen ce between our results and those of 
IGilbank et al.l (|2008l ) seem to be as they go in opposite di- 
rections. Unfortunately it is difficult to m ake a straight- 
forward comparison between our results as IGilbank et al] 
((2008) split their clusters at a richness that corresponds 
to ~ 400 km/s, below which we only have four clusters. 
We therefore have almost no constraint on the behavior 
of iVium/-^faint for thcsc vcry poor systems. To resolve 
this discrepancy will require computing iViumZ-^faint in 
the Gilbank sampl e with our velocity dis persion cut. An- 
other difference in IGilbank et al.l 1)20 08') with respect to 
our work is that those authors use a color cut that iso- 
lates only the red side of the red sequence, whereas our 
color cut encompasses a band centered on the peak of the 
red sequence. It is possible, therefore, that our cut has a 
different contribution from blue galaxies at the blue side 
of the red sequence. It is not clear, howe ver, which cut 
is mor e physically meaning ful. While the IGilbank et all 
(|2008l ) cut only isolates galaxies with the reddest colors 
it may miss galaxies that have been most recently added 
to the red sequence and therefore those that have the 
youngest stellar population ages and bluest colors. At 
present we cannot determine the nature of the discrep- 
ancy as a cut only on the red side of the red sequence 
in our data would substantially reduce the significance 
which we can measure A^ium/-^faint- 

We do not have enough galaxies to reliably fit a and 
M* independently for our composite LFs split by veloc- 
ity dispersion. Because these two parameters are highly 
degenerate and there are indications that Mum/A^faint is 
different in the different bins of cr we therefore do not 
attempt to interpret the Schechter function fits split by 
a. 

5.2.3. radial gradients 

We examined whether there was any evidence for radial 
gradients in the LF. We do not have enough galaxies to 
split our clusters up into annuli so instead we compared 
the LFs within 0.75 and 0.5 Mpc, with the acknowledg- 
ment that these two are correlated. For both the clusters 
at 0.4 < z < 0.6 and those at 0.6 < z < 0.8 the galax- 
ies at 0.5 Mpc < r < 0.75 Mpc make up 30% of the 
galaxies at r < 0.75 Mpc. For the EDisCS clusters we 
made a similar comparison as in Figure [Hj and find that 
there are no significant differences in the LFs at different 
radii. As expected from this, there are also no significant 
differences in the Schechter function parameters. These 
results are tr ue for all rest-fr ame filters. In their study of 
local clusters. fPaolillo et al.l (,2001) find results consistent 
with ours but for the whole galaxy population ins t ead o f 
just for those on the red sequence. iPopesso et al.l ()2006D 
measure the red sequence LF at different radii in local 
clusters and find that there are radial trends in the LF, 
but primarily at magnitudes fainter than what we probe 
with our data. At magnitudes corresponding to those 



we probe they find only a weak dependence of the LF 
with cluster radius, which is entirely compatible with 
our lower precision measuremen t of the radial trends. 

1 Hansen et al.l (jsubmitted to Apjl ) measure the LF of red 
galax ies in different radia l bins in many more clusters 
than iPopesso et al.l ()2006f ) but to significantly brighter 
magnitude limits and also find no radial dependence in 
the shape of the LF. 

6. DISCUSSION 

6.1. Mass dependent build-up of red sequence cluster 
galaxies 

As shown in i i5.2.1l and Figure |5| there is dramatic evo- 
lution in the LFs of red sequence galaxies from z = 0.8 
to z = 0. At the same time this evolution appears to 
proceed at different rates for galaxies of different lumi- 
nosities. Red sequence galaxies at L > L* appear to 
be in place at z ~ 0.8 but the fainter galaxies were 
added to the red sequence at much later times. As a 
different way to visualize this we plot in Figure |8| the 
ratio of the passive-evolution-corrccted EDisCS LFs to 
those from the SDSS. At bright magnitudes both LFs 
are consistent with a constant ratio, implying that the 
EDisCS and SDSS LFs have the same shape. At fainter 
magnitudes the ratio then decreases. In the intermedi- 
ate redshift clusters it appears that there is a threshold 
magnitude, fainter than which there is a deficit of red 
sequence galaxies with respect to the SDSS clusters and 
above which there appears to be a constant ratio. For 
the high redshift clusters there is not enough signal-to- 
noise to determine if a threshold magnitude exists or if 
there is simply a monotonic trend to fainter magnitudes. 
In any case it is clear that the magnitude brighter than 
which the red sequence is in place is fainter at lower red- 
shifts. This is reminiscent of the results bv iBundv et al.l 
(j2006) who show that there is a threshold stellar mass 
above which star formation appears to be quenched in 
field galaxies, and that this threshold evolves to lower 
masses at lower redshifts. However, cluster-associated 
processes for quenching star formation, e.g. ram-pressure 
stripping, strangulation/starvation, or harassment, are 
different from the mass-dependent quenching that may 
be present in the field. It may be that the evolving mass 
threshold in the field is imprinted on the cluster red se- 
quence in the form of infalling red galaxies whose star 
formation is halted for reasons not associated with the 
cluster. We cannot, however, directly confirm that hy- 
pothesis with our data. It will be very interesting in the 
future to compare the mass functions of cluster and field 
galaxies to see if there are differences between the cluster 
and field. 

These results imply that the bright cluster galax- 
ies were already in place and on the red sequence by 
z = 0.8. This is consistent with the evolution of the 
colors of these galaxies, which is fit well by a SSP with 

2 < Zform < 3 (Standford et al. 1998; Holden et al. 
2004; DL07). It is also perfectly consistent with the 
formation redshifts found by fundamental plane stud- 
ies (e.g. van Dokkum & Franx 2001; van der Wei et 
al. 2005; J0rgensen et al. 2006; van Dokk um fc van der 
Mare] 2007; Saglia et al. in preparation). iBarger et al.l 
(jl998) found no formal evidence of evolution in their K- 
band LF for morphologically classified early type galax- 
ies, but had large enough error bars to be compati- 
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Fig. 6. — The ratio of the LF of red sequence galaxies in clusters of different velocity dispersion. Large differences in the shapes of 
the LFs will appear as significant deviations from a horizontal line in this plot. Each column indicates a different redshift bin and each 
row represents a different rest-frame band. The ratio has been normalized to a mean of unity to demonstrate relative differences with 
magnitude. For the g-band LF for EDisCS clusters we show the ratio of LFgs as red triangles to compare to the ratio of LFzp given as black 
squares. All EDisCS points for the r and j-bands are usi ng LFzp. The v e rtical dotted lines in the bottom row give the division between 
bright and faint galaxies that correspond to those used in IDe Lucia et al.l I I2007I) . 

have shown that there is indeed a differential build-up of 
the red sequence with redshift, such that fainter galaxies 
were added to the red sequence at lower redshifts than 
brighter galaxies. This may in fact be consistent with 
claims that the early type fraction increased significantly 
at z < 1 (e.g. Dressier et al. 1997; Fasano et al. 2000; 
Smith et al. 2005; Postman et al. 2005; Desai et al. 
07) implying a late buildup of significant fractions of the 
red galaxy populat ion. Th is is explored in more detail 
in Sanchcz-Blazqu ez et al.l ([2009), which examines the 
joint evolution of ages, metallicities, and morphologies 
of galaxies on the red sequence. 

6.2. Comparison to field 

We compare our red-sequence LF to two independently 
computed red sequence LFs for the field^** . The first field 
LF is that pubhshed bv lBrown et all ()2007D in the NO AO 

^* There are other measures of the early-type galaxy LF in the 
field at z = f (e.g. Zucca et al. 2006; Scarlata et al. 2007) but 
we choose not to com pare directly with these surveys for various 
reasons. For example, IZucca et al.l II2006I ) identify red galaxies by 



ble with the expected amount of passive evolution over 
their redshift interval. They did perform an analysis 
of the surface brightness evolution, however, and found 
the amount of fading expected from passive evolution. 
Later works per formed more extensive LF analyses us- 
ing K-hand LF (Ide Propris e t al. 19 91), and i n the rest- 
frame NIR (|Andreonll2006t iDe Propris et al.l [20071 ) and 
find evolution in M* consistent with passive evolution 
and Zforra > 1-5 when fit to galaxies with a faint limit 
1 — 3 magnitudes brighter than ours. This may imply 
that the red sequence galaxies dominate the galaxy mass 
function of clusters even out to z ~ 0.8, such that the 
total LF of bright galaxies in the rest-frame NIR is re- 
ally dominated by red sequence galaxies. We will address 
this directly in our future analysis of the galaxy stellar 
mass function in our clusters (Aragon-Salamanca et al. 
in preparation). 

In previous works, the high inferred formation redshift 
for bright red galaxies has been used by many LF studies 
to imply that the population of all red sequence cluster 
galaxies was already in place at 2; > 1. We, however. 
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Fig. 8. — The ratio of the EDisCS composite red sequence LFs to 
those from the SDSS in the rest-frame r, and j-bands. In all cases 
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intermediate and high redshift EDisCS clusters. The horizontal 
lines is drawn at unity to guide the eye. In all three bands the 
bright end is consistent with a constant value, which appears to 
turn over at fainter magnitudes. The turnover magnitude evolves 
to fainter magnitudes at lower redshifts. 



spectroscopic type and it is not clear how that corresponds to our 
selection by color. IScarlata et al.l II2007I) do not aperture correct 
their magnitude measures and their LFs are 1-2 magnitudes shal- 
lower than ours, while not having a higher signal-to-noisc than the 
NDWFS. 



Deep Wide Field Survey (NDWFS). The NDWFS LF is 
computed over 6.96 deg^. It is somewhat shallower than 
the EDisCS survey but the LF has excellent signal-to- 
noise at both the bright and faint-end of the LF. Magni- 
tudes in the NDWFS were measured in circular apertures 
of r = 4'.'0 and had an aperture correction derived from 
simulations. Although our surveys have substantially dif- 
ferent depth and image quality, we have checked that our 
total magnitude estimates are compatible with those in 
the NDWFS. We compared the aperture corrections we 
apply to our AUTO magnitudes in EDisCS to the differ- 
ence between the AUTO magnitudes of NDWFS galaxie s 
and their total magnitude estimate (jBrown et al.ll2007b . 
This difference is similar, indicating that our correction 
brings our total magnitudes into rough agreement with 
theirs. Both our survey and the NDWFS select red galax- 
ies in similar ways, implying that our results can be di- 
rectly compared. 

For our second comparison we determine the LF for 
red sequence galaxies using a modified ver sion of t he data 
and LF technique presented in [Ma rche sini et al.l ([2007). 
The original data are comprised of multi-band optical 
through NIR photometry over 6 dispersed fields, with a 
total area of 355 arcmin^. The fields are comprised of 
the Hubble Deep Field South, The Chandra Deep Field 
South, and the four fields from the Deep NIR Multi- 
Wavelength Survey by Yale-Chile (MUSYC; Quadri et 
al. 2007). The main modification that we ha ve made to 
the data presented in lMarchesini et all (j2007[ l is that we 
have included photometry past the ii'-band to 8.0/im us- 
ing Spitzer observations (Wuyts et al. 2008; Marchesini 
et al. 2009). We have also determined the photometric 
redshifts using a different code (EAZY; Brammer et al. 
2008). The data are if-selected but we used the observed 
I — K colors of galaxies at z < 1 in the EDisCS /-band 
selected catalog to verify that we can detect all red se- 
quence galaxies in the K-hand data at z < 1. There- 
fore we can use the isT-band selected data to construct a 
pseudo /-band selected red sequence sample analogous to 
that for EDisCS. Total magnitudes in the MUSYC sur- 
vey are constructed using SExtractor AUTO apertures 
with aperture corrections derived in an identical way as 
our own. The MUSYC red sequ ence definition is id en- 
tical to our own and that from iBrown et al.l (|2007l ) is 
very similar, ensuring that we are selecting red sequence 
galaxies in the same way. 

The NDWFS red gal axy LFs are co nsistent with 
those from COMBO-17 (iBell et al.l [200l and DEEP2 
(|Willmer et "all [20061 : iFaber et al.ll200l " but are derived 
over a much larger area and have a very well deter- 
mined bright end. The NDWFS LF is only computed 
in the rest-frame S-band and so we compare it to the 
EDisCS LF also computed in the rest-frame i?-band. 
The MUSYC LFs are statistically consistent with the 
NDWFS LFs at bright magnitudes but go significantly 
deeper and have LFs in the the g, r, and i-bands. In both 
field samples, the Schechter function fits to the field data 
have been performed using a maximum-likelihood tech- 
nique with a as a free parameter and the plotted points 
are the 1/T4nax measurements of the LF which include the 
contribution from field-to-field variance for the MUSYC 
data. 

We compare the field LFs to those from EDisCS in 
Figure [9] and [101 In both figures for display purposes 
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we have normalized the different LFs to have the same 
integrated luminosity. The bright ends of the field and 
clusters are consistent in all cases. At fainter magnitudes 
the cluster LF appears be be over abundant compared to 
the field at 0.4 < z < 0.6 but this trend reverses itself 
at 0.6 < z < 0.8. The reversal is most apparent at the 
reddest rest-frame wavelengths (the right two columns 
of Figure [To)) , where the MUSYC LF extends to similar 
depths as the EDisCS LF. 

With regards to the best-fit Schechter function param- 
eters, those from the two field surveys agree with each 
other at better than 1.5-sigma. It appears that the field 
value of M* is systematically brighter than for the clus- 
ters in both redshift bins, but only at the 2-sigma 
level. At 0.4 < z < 0.6 the faint end slope of the clus- 
ters is steeper than for the field at the ^ 2-sigma level. 
At 0.6 < 2 < 0.8 however, this has reversed, with the 
field having an a that is ~ 2 sigma steeper than for 
clusters. The change in the relative slopes is interesting 
as, in addition to the dependence of Mum/A'faint evolu- 
tion on cluster velocity dispersion presented in ij5.2.21 it 
provides further indication that the buildup of the red 
sequence happens at different rates in different environ- 
ments. Further it appears that the population of the 
faint-end of the red sequence happens more quickly in 
clusters than in the field. This could occur if the main 
channels for the population of the red sequence in the 
field occur at a roughly constant or decreasing rate with 
increasing cosmic time, as may be the case for AGN feed- 
back (e.g. Croton et al. 2006) or galaxy-galaxy mergers 
(e.g. Lotz et al. 2008), while quenching processes associ- 
ated with clusters become more efficient with increasing 
time. This is plausible in a ACDM universe where cluster 
assembly proceeds most rapidly at late times, implying 
that the processes associated with quenching in clusters 
would also become more efficient as one moves towards 
lower red shift. At face valu e, this is consistent with the 
results ofl Pesai et al.l ()2007t) . who fo und that the SO frac - 
tion in clusters from EDisGS and iFasano et al.l ()2000r ) 
was changing slowly at z > 0.5 but increased much more 
rapidly starting dX z — 0.4 — 0.5. While, the exact time 
at which the red sequence assembly in clusters and the 
field cross and the time at which the SO buildup be- 
comes significant are not very well determined it is pos- 
sible that the difference in redshift is real. If it is we 
may speculate that the 0.7 Gyr difference in time be- 
tween z — 0.6 and z — 0.5 may reflect an intrinsic delay 
between morphological transformation and the trunca- 
tion of the SFR. As an example of such a scenario, the 
SFR may have been truncated during one pass through 
the cluster, while the morphological transformation may 
have r equired repeated cluster pass ages to build up the 
bulge (iGhristlein fc Zabludofil[200l. 

It is worth noting that iGilbank fc Baloghl (|2008f l com- 
pile many different measurements of Mum/A^faint in the 
field and cluster and find that Aium/A'^faint evolves more 
quickly in the field than in the cluster. It is difficult to 
directly compare our results as we do not directly com- 
pare Mum/A^faint between the field and cluster. We also 
find that the difference between the field and cluster is 
most dramatic at red re st-frame wavelengths, whereas 
IGilbank fc Baloghl (12008") measured their Aium/A^faint in 
the rest-frame T^-band. It is also curiou s that the field 
Mum/A^faint in IGilbank fc Baloghl (|2008l ) is observed to 



evolve quicker than for clusters but in DL07 and in this 
work, high-mass clusters are observed to evolve quicker 
in Aium/A'^faint than low-mass clusters (see ij5.2.2l for a 
discussion of these results). All we can say with some 
certainty is that the bright ends of the field and cluster 
measured LFs agree within our, admittedly, large error- 
bars but that the faint ends do not, that this disagree- 
ment increases towards redder rest-frame bands, and that 
there seems to be a tentative indication that the direction 
of this disagreement changes over the EDisGS redshift 
range. 

It is also interesting to discuss our results in light of 
similar cluster vs. field comparisons in the local Uni- 
verse. Our finding that the faint-end slope in clus- 
ters is steeper than in the field at 0.4 < z < 0.6 
is in qualitative agreement with the local results using 
spect roscopically defined no n-starformi ng galaxies from 
2dF (de~Propris et al.ll2003[) . However. Ide Propris et al.l 
(12003 ) found that the clusters have a brighter M* than 
the field, which may be inconsistent with our results at 
0.4 < z < 0.6. This may imply the presence of rela- 
tive evolution in the bright end of the LF between the 
field and the cluster but it is important to keep in mind 
that our composite cluster LF is noisy at the bright 
end and that the local studies used spectroscopic tech- 
niques to identify galaxies with no star formation. Keep- 
ing this in mind, as we showed in Figure [5] the bright 
end of the measured EDisGS LFs are consistent with 
pure passive evolution, which would argue against a sig- 
nificant increase in the cluster red sequence population 
at the luminous end toward lower redshift. One pos- 
sible explanation is that red galaxies in the field are 
younger than those in clusters, and will therefore fade 
by a larger amount toward lower redshift. A useful 
check of this c omes from Fundamental Plane studies. 
Ivan der Wei et a l. (2005) have shown that the evolution 
in M/Lb for massive galaxies in the cluster and field is 
Aln{M/LB) = (-1.12 ± 0.06)z and (-1.2 ± 0.18)z re- 
spectively. These correspond to AMb = (— 1.22 ± 0.06)z 
and (— 1.30 ± 0.20)z for cluster and field galaxies respec- 
tively. The lack of a difference between the LF evolution 
of red galaxies in the field and in c lusters agrees with 
Ivan Dokkum fc van der Marell ()2007f ) who have shown 
that massive ellipticals in clusters and the field have lu- 
minosity weighted ages t hat are within 4.1% (ss 0.4 Gyr) 
of each other. Gcbhar dt et al.l ()2003D . however, find a 
~ 2 Gyr difference in th e age of cluster and field ellip- 
ticals but the analysis of Ivan Dokkum fc van der Marell 
(|2007f ) involves more sophisticated dynamical modeling 
and also co rrects for selection effects su ch as, "progen- 
itor bias" (Ivan Dokkum fc Franxl 1200 ll ). In any case, 
this small difference between clusters and the field is not 
en ough to explain the diffe rence that we see with respect 
to Ide Propris et al.l (|2003f l. Another explanation must 
therefore be found to explain the apparent difference be- 
tween the M* in the field and clusters at intermediate 
and low redshifts and what it implies for the relative 
evolution of bright ellipticals in these two extremes of 
environment. To better study the relative evolution of 
the bright end of the LF in the field and in the cluster it 
will be necessary to construct red sequence LFs for much 
larger samples of clusters. 

6.3. The integrated growth of the red sequence 
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Fig. 9. — A comparison of the B-band LFs and Schechter func- 
tion fits for the composite LP of red sequence galaxies in EDisCS 
clusters and the LF of red sequence galaxies in the field as deter- 
mined from the NDWFS (Brown et al. 2007) and from MUSYC, in 
two redshift bins. The vertical arrows give the values of M* and its 
68% confidence limits. The upper, middle, and lower arrow refer 
to the MUSYC, EDisCS, and NDWFS LFs respectively. We have 
used the NDWFS Schechter fits that allow a to vary. The field 
LFs have been scaled vertically to have the same total luminosity 
density as the EDisCS LF. 
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Fig. 10. — A comparison of the g, r, and i-band LFs and 
Schechter function fits for the composite LF of red sequence galax- 
ies in EDisCS clusters and the LF of red sequence galaxies in the 
field (as determined from the MUSYC survey), in two redshift 
bins. The vertical arrows give the values of M* and its 68% con- 
fidence limits. The upper and lower arrow refer to the EDisCS 
and MUSYC survey respectively. The field LF has been scaled 
vertically to have the same total luminosity density as the EDisCS 
LF. 



In this subsection we explore the buildup of the total 
amount of light on the red sequence as the clusters evolve 
from z ~ 0.4 — 0.8 to z ~ 0. We start by measuring the 
total light on the cluster red sequence, jcis, by integrat- 
ing the measured LFs in each cluster. We do this both 
at r < 0.75 Mpc and at r < 0.5i?2oo: which scales with 
M200- We present the results here for r < 0.75 Mpc. The 
choice of aperture for the jcrs computation is somewhat 
arbitrary, but we note that the results using r < 0.5i?2oo 
are consistent with those computed using r < 0.75 Mpc 
but the trends are not as significant. In the left-hand col- 
umn of Figure [Tl] we compare the total light on the red 
sequence in each EDisCS cluster calculated both using 
the measured LF down to the magnitude limit of each 
cluster and using the complete integral of the best-fit 
Schechter function. The right-hand column is the same 
for blue galaxies, which we discuss later in this section. 
There are 3-4 clusters (with the number depending on 
the rest-frame bandpass) which have low measured lu- 
minosity and as a consequence also did not have well 
constrained Schechter fits (e.g. CL1227. 9-1138 in Fig- 
ure |4]). For the rest of the clusters, however, the two 
measures of jcrs correlate very highly, with a mean offset 
of 3 — 5% and a scatter of 9%. This demonstrates that 
our observations go deep enough to directly probe almost 
all of the light on the red sequence. Although the formal 
errors on the integral of the Schechter functions are much 
smaller than for the measured LF integrals, this is pri- 
marily because of our assumption of a parametric form 
for the LF and because we assume a well constrained 
faint-end slope, as determined by a fit to the composite 
LF in each redshift bin. To be as conservative as possible 
we therefore use jcrs derived from the measured LFs for 
the rest of this discussion. This has the added advantage 
of allowing us to include some of the lowest luminosity 
clusters that had poorly constrained Schechter function 
fits. 

To perform a consistent comparison between low and 
high redshift clusters, and to mitigate any secondary de- 
pendences on cluster mass we compared our clusters to 
those in the local universe as a function of cluster mass. 
We calcu late the mass, M2n o, from the velocity dispersion 
following [Finn et all (|2005( ): 
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In Figure [T^ we plot jcrs vs. M200 for the EDisCS clus- 
ters and for the SDSS clusters. Even at a fixed mass in 
the SDSS sample there is significant scatter in integrated 
luminosity. This is intrinsic scatter in the cluster pop- 
ulation as ^ 95% of the SDSS clusters have lower than 
20% error on the integrated luminosity. Although we 
have a large sample of SDSS clusters, the size of the lo- 
cal sample, especially at the massive end is the dominant 
uncertainty in the following analysis. The EDisCS clus- 
ters lie at brighter luminosities than the SDSS clusters at 
the same cluster mass. At a basic level this is expected 
because the red galaxies will fade as a function of time, 
moving the EDisCS points in the direction of the SDSS 
locus. We will quantify this evolution below. 

It is also clear that the SDSS clusters deviate signif- 
icantly from the one-to-one relation between M200 and 



20 



Red Sequence Cluster Galaxy LF at z < 0.8 



jcrs(i-e- constant jcrs/M2oo), which indicates that there 
is a residual dependence of jcrs/M2oo on M2oo- In this 
case it appears that more massive clusters have a lower 
jcrs/M2oo than less massive clusters. We have confirmed 
that this trend is not due to the uncertainties on the ve- 
locity dispersion for SDSS clusters. A similar deviation 
from the one-to-one relation between jcrs and M200 ex- 
ists for the EDisCS clusters, such that clusters with low 
M200 are brighter than the one-to-one relation, but it is 
not clear how robust this is given the small numbers of 
very low mass clusters. 

In Figure [Ol we plot jcrs/M2oo vs. M200 for the EDisCS 
and SDSS clusters. In the left column we plot the individ- 
ual values and in the right column we show the geometric 
mean jcrs/M2oo of the EDisCS clusters in three mass bins 
and compare them to the SDSS clusters. The mass bins 
were chosen to contain approximately similar numbers 
of objects. We do not have enough clusters to bin both 
in mass and in redshift, but there is no dependence of 
sigma and redshift in our sample (see Figure I14p and 
a Spearman rank test gives only a 4% probability that 
the two values are correlated, i.e. our mass bins should 
not be affected by secondary correlation of cluster mass 
with redshift. Again, as in Figure [T^ it is clear that the 
jcrs/M2oo values for the EDisCS clusters are larger than 
for the SDSS clusters. 

Even though the EDisCS clusters have higher jcrs for 
their M200 s.t z 0.6 they must by definition evolve 
by 2; onto the local relation defined by the SDSS 
clusters. To explore what this constraint implies for 
the build-up of light (or mass) on the red sequence we 
have constructed a set of 4 toy models which we de- 
scribe below. In all models we assume that two processes 
are universally present. First, the mass of clusters will 
grow via accretion of matter from the surrounding cos- 
mic web. Second, the galaxies on the red sequence at 
the epoch of observation will only fade as a function of 
time, i.e. we assume that galaxies on the red sequence 
are passive and stay on the red sequence. We account 
for the growth of clusters by tracking the median growth 
in mass for cluste r s of a cert ain mass using the r esults 
of iWechsler eTall (l2Q02h and iBullock eTallllMlh ^^'^s. 
To estimate the fading of cluster red sequence galaxies 
we fit a linear relation to the change of magnitude as a 
function of log(time) for each bandpass. In Figure [T5] 
we show the fading in magnitude for models of two dif- 
ferent metallicities ()Bruzual fc CharlotI [2003f ). A linear 
fading of magnitude with log(age) is a very good ap- 
proximation for SSPs with age> 1.4Gyr, as was shown 
originally bv iTinslevI ()1980l ). Models with Z=Zq and 
2.5Zo differ in the amount of fading by < 0.13 magni- 
tudes over the 5.7Gyr period from z = 0.6 to 2 = 0. We 
adopt the Tj—Tjq model but our results are insensitive to 
this choice. In tracking the fading stellar populations we 
evaluate whether galaxies fade below the absolute mag- 
nitude limit that we adopt for the SDSS. The first, sim- 
plest, and most unrealistic toy model is one in which 
only mass is accreted but no light is added to the red 

We computed this using programs obtained from 
http : //www.physlcs .uci . edu/~bullock/CVIR/ 

Note that our highest mass clusters will evolve into objects 
that are so rare as not to be present at all in the SDSS C4 sample. 
This is also the fate of most of the numerous massive X-ray clusters 
at 2 > 0.5 that are found in the literature. 



sequence of the cluster. The light in the cluster there- 
fore decreases by pure fading. In the right-hand panel of 
Figure [T^ we demonstrate how this model would cause 
the EDisCS clusters to evolve from Zdust down to z = 0. 
This model results in predicted jcrs/M2oo values at z = 
that are too low compared to the SDSS. In Figure [16] 
we demonstrate this in another way by plotting the ratio 
of the predicted z = 5-band red sequence luminos- 
ity, jcrs,pred^ to the observed SDSS red sequence lumi- 
nosity jcrs,SDSS- From this figure it is clear that this 
first model results in predicted z = luminosities that 
are too faint by a factor of 1 — 3 depending on the 
mass range, although formally the required mass growth 
implied in all three mass bins are consistent with each 
other. The red sequence light in clusters must therefore 
grow by a similar factor from z ~ 0.6 to the z ~ 0. For 
clusters with M2oo< lO^'^'^MQthe required growth is a 
factor of 2.5-3. Because light and mass are proportional 
for red sequence galaxies, this also implies that the stel- 
lar mass on the red sequence needs to therefore grow 
by a factor of 1 — 3 with a growth of 2.5-3 required at 
M2oo< 1O^^'^A^0. That the largest implied growth may 
come from clusters of low to intermediate masses rn ay be 
consistent with the results of lPoggianti et al.l ()2006[ ). who 
find that the star-forming fraction in clusters at z < 0.8 
evolves most rapidly at intermediate cluster velocity dis- 
persions (M200) but evolves little at the highest velocity 
dispersions. 

To explore different scenarios for how additional light 
(and stellar mass) may be added to the red sequence 
we therefore consider a second model in which the clus- 
ter mass growth is accompanied by the accretion of 
stars with the same Mtot/L and its expected evolution, 
where Mtot/L is derived for each cluster individually 
using the total light in red sequence cluster galaxies. 
In other words, this conservative model assumes that 
the only galaxies added to the red sequence are red 
sequence field galaxies with the same ages and SFHs 
as cluster ellipticals and with the same proportion of 
light relative to the total mass. This model would be 
consistent, for example, with the very small age differ- 
en ces (~ 4%) found between clus t er and field ellipticals 
by Ivan Dokkum fc van der Marell ()2007l ) . Assuming all 
galaxies that fall in to clusters at z < 0.6 end up on the 
red sequence at z = this model represents the minimum 
amount of light that can be added to the cluster red se- 
quence by accretion from the field. The results of this 
model are shown as the solid triangles in Figure [TH As 
expected these models yield higher predicted z = lumi- 
nosities for the EDisCS clusters, and are more consistent 
with the expectations for SDSS clusters. 

It is interesting to discuss the predictions of this model 
in relation to the roughly factor of two grow th in mass 
on the red sequence inferr e d from field s t udiesjBell et all 
(|2004[ ): lFaber et al.l (|2007l ): lBrown et al](|2007l ). In inter- 
preting this it is important to remember that the "field" 
surveys contain a range of environments, including mod- 
erate mass clusters for the largest surveys like the ND- 
WFS and extending to massive groups for the MUSYC, 
DEEP- 2, and COMB 0-1 7 surveys. The observed mass 
growth on the red sequence in field surveys represents 
the actual transformation of blue galaxies to red galax- 
ies. Clusters, on the other hand, are growing their to- 
tal mass by a factor of '--^ 2 at z < 0.6 and, at least 
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partly, will be increasing the total amount of red light 
merely via the accretion of red galaxies from the field 
as in model 2 (above). Thus, the total increase that we 
infer in the mass on the red sequence is consistent with 
only a moderate additional transformation of blue to red 
galaxies. Nonetheless, our value for the required mass 
growth is still rather uncertain due to the significant 
intrinsic scatter both in the SDSS and EDisCS cluster 
jcis/M2oo values. Also, our middle and lowest mass bins 
(M2oo< imply a factor of --2.5-3 growth 

in the red sequence stellar mass at z < 0.8, which may 
exceed the observed growth in the field and the expected 
total mass growth in clusters. As we show in §6.2l there is 
differential evolution in the shapes of the LFs in the field 
and clusters and this does imply that clusters evolve more 
rapidly in the number of g alaxies on the red sequ ence. 
There is also evidence from iPoggianti et al.l (|2006f ) that 
star forming galaxies are be ing truncated pre ferentially 
in cluster environments but iFinn et aTl ()2008[ ) find that 
the decline in the SFRs of cluster galaxies at 2 < 0.8 is 
comparable to that in the field. While it is difficult to 
draw precise conclusions about the necessary mass trans- 
fer to the red sequence, it is clear that a pure passive 
fading of the cluster red sequence seen at z = 0.4 — 0.8 
will result in clusters that are systematically too faint 
compared to those seen locally. Within the large uncer- 
tainties in the model predictions there is no significant 
trend with cluster mass, which is also true for each of the 
following models. 

As a third model we calculate how much light is added 
to the EDisCS cluster red sequences by z = by all of 
the blue galaxies in the clusters at z ~ 0.4 — 0.8, with no 
additional infall of red galaxies (pentagons; Figure [TB]) . 
In calculating this estimate we take into account the un- 
certainties resulting from the differences between LF^p 
and LFss that we determined in §4.51 Indeed, the total 
luminosity of blue galaxies in the g-band, jbg, determined 
from LFss ranges from 0.7 to 5.7 times larger than jtg de- 
termined from LFzp, with a median of 1.8, as calculated 
over all clusters. In addition, as is shown in the right- 
hand column of Figure [Tl] the blue galaxy g-band LFgs 
has not converged completely for low luminosity clusters 
due to the very steep faint-end slope. Nonetheless, for 
most of the clusters the missing light below our magni- 
tude limit is small and jbg as derived from LFgs and LF^p 
should still bracket the true value of jbg ■ For every clus- 
ter we therefore use the mean of jbg as determined from 
LFss and LF^p and the values for the two methods as an 
indication of the uncertainty in jbg ■ Because we can only 
calculate both LFgs and LFzp for the g-band (see 
we limit our analyses for the following models to that sin- 
gle bandpass. We then assume that the galaxies in each 
cluster have been forming stars constantly prior to the 
epoch of observation and that they have achieved solar 
metallicity by the time they are observed. We assume 
that they have Ay — 1, which corresponds to Ag = 1.17. 
Subsequent to observation we assume that these galaxies 
continue forming stars for 1 Gyr before abruptly ceasing 
their star formation and losing their dust. Except for the 
extra extinction, this is similar to the delayed truncation 
model of DL07 and is consistent with the evolution in the 
shape of the LF from z = 0.8 to 0. For every cluster we 
calculate the expected luminosity contribution that these 
blue galaxies make to jcrs,pred assuming that the galaxies 



were forming stars for 3~6 Gyr prior to observation. As 
we can see in Figure [16] this simple model over predicts 
the amount of light on the red sequence in local clusters, 
especially so for the highest and lowest mass bins. 

In reality, galaxies must be accreted onto the cluster 
over time and we therefore consider a fourth model in 
which both old ellipticals are accreted onto the cluster 
(model 2) and blue galaxies within the cluster are trans- 
formed (model 3). Being the sum of models 2 and 3, 
this fourth model naturally also over predicts the pre- 
dicted light in local clusters, by a factor of 1.8-3.6. In 
understanding why this last model over predicts the local 
luminosity of clusters it is important first to remember 
that these models are in many ways very conservative in 
how much light is added to the red sequence by z = 0. 
The amount of mass that our clusters accrete is deter- 
mined by ACDM and for this accreted mass we add the 
smallest possible amount of light to the red sequence by 
only accreting old galaxies at the same jcrs/M2oo as the 
cluster. Accreting galaxies that are not as old as cluster 
ellipticals (e.g. blue galaxies) will increase the luminosi- 
ties of the clusters by z = with respect to our fourth 
model. In addition, there is little way to avoid the trans- 
formation of blue cluster galaxies at z = 0.6 to red and 
dead ones by z = so our third model should also be 
valid. Because we accrete as little light as possible for 
the expected mass accretion, our fourth model can be 
thought of a lower limit on the amount of light added 
to the cluster (but see below). In light of these rather 
conservative assumptions it is perhaps puzzling that the 
simple model nonetheless produces too much cluster red 
sequence light at z = 0. As a note, the amount by which 
model 3 and 4 over predict local luminosities of clusters 
when computing jcrs at r < 0.5i?2oo is consistent with, 
but not as large as when computing jcrs at r < 0.75 Mpc. 

There are at least three possible resolutions to this 
apparent discrepancy. First, it may be that clusters ac- 
creted significant light in galaxies that never enter the 
red sequence by z = 0, i.e. from star-forming galaxies 
that have not had time to cease star formation and red- 
den since they fell in. In our SDSS sample we also com- 
puted the blue luminosity function and find that 26 ±3% 
of the total cluster light comes from galaxies bluer than 
the red sequence. Reducing the contribution to the fi- 
nal red sequence luminosity by this amount corrects for 
our assumption that all accreted galaxies end up on the 
red sequence. This brings the fourth model into closer 
agreement with local cluster luminosities but still sys- 
tematically over predicts them. 

A second possibility is that blue star-forming galaxies 
at z = 0.6 that enter the red sequence at z = may still 
be dust enshrouded, lowering their total observed lumi- 
nosity, som ething not encapsu lated in our simple mod- 
els. Indeed IWoIfeTall (|2005h find that roughly 1/3 of 
red sequence galaxies in the Abell 901/902 superclus- 
ter are dusty star forming galaxies with (Ay) « 0.6. If 
the same fraction of dusty red galaxies is present in all 
clusters, and under the simplifying assumption that the 
fraction of dusty galaxies is independent of galaxy lu- 
minosity, then this extinction changes the jds by a fac- 
tor of 0.87. Taken alone this is obviously too small to 
make a significant contribution to reconciling model 4 
with the local jcrs values. However, combining this with 
the amount of light from accreted galaxies that never 
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make it onto the red sequence (see previous paragraph) 
would change jcrs,pred by a factor of 0.64, which still re- 
sults in jcrs,pred valucs that are systematically a factor 
of 1.15-2.3 too high compared to local clusters, but are 
consistent within the 68% confidence limits. 

Thirdly, it is possible that a substantial fraction of 
stars / mass is in cluster components other than red 
sequence galaxies, and is thus neglected in our mea- 
surement of the LF. For example, our LFs exclude 
the BCG and the intracluster light (ICL) which could 
harbor a significant fraction of the total stellar mass 
and may evolve differently than the red galaxies. In- 
deed, cluster galaxies are expected to sink to the clus- 
ter center and merge wi t h the BCG ( "cannibalism" 
lOstriker fc Tremai nelll975t IWhitd [19761 ). a phenomenon 
that has been observed in some low redshift BCGs fiLauer. 
119881 ). Cannibahsm may be especially relevant for the 
most massive cluster galaxies, for which the dynamical 
friction timescale is similar to the timescale we are prob- 
ing here (a few Gyrs). Hence, the extra mass in stars our 
toy model predicts c ould simply have bee n accreted onto 
the BCGs. However. IWhilev et all (|2008l ) found that the 
properties of the EDisCS and SDSS BCGs are consistent 
with passive evolution since z 2, implying no apprecia- 
ble BCG mass growth. This result seems at odds with 
the fact or 3-4 mass growth of BCGs predicted by simu- 
lations (|De Lucia fc Blaizot 20071). One needs to keep in 
mind, however, that .Whilev et all (|2008D considered only 
the central 37 kpc of each BCG and that any mass ac- 
creted in mergers must predominantly be accreted onto 
the envelope of the BCG and/or the ICL. This is also a 
viable possibility to explain the discrepancy between our 
toy model prediction for the mass on the red sequence: 
rather than remaining in the galaxies we see in the clus- 
ter at z ~ 0.6, a significant fraction of old stars could 
be part of the ICL at z = 0, and are thus not accounted 
for in our LFs. There are various mechanisms by which 
cluster galaxies may get disrupted and lose their stars to 
the ICL: e.g. strippin g due to tidal forces in the clus- 
ter gravitational field (Merritt '1984), or galaxy harass- 
ment (Farouki & Shapiro 1981; Moore et al. 1996). In- 
deed, the mass in the ICL has been measured i n a nu mber 
of low-redshift clusters, e.g. iGonzalez et al.l ()2007f ) find 
that the BCG+ICL contribute 20%-40% of the stellar 
light within i?5oo in clusters of the range in velocity dis- 
persion we are co nsidering her e . Th is result is consis- 
tent with that of IZibetti et al.l (|2005l ). who found that 
the BCG-I-ICL contribute ~30% of the stellar light in 
stacked SDSS clusters. Furthermore, the color of the 
ICL is comparable to, or ev en slightly redder t han, the 
total color of cluster galaxies lZibetti et al.l ()2005f ). imply- 
ing that t he ICL may hav e originated from red sequence 
galaxies. IGonzalez et al.l ((200 5i) use the position angle 
and ellipticity for a set of lower redshift clusters to de- 
compose the BCG and ICL and find that 80% of the 
light on average comes from the ICL. We do not know 
the ICL contribution in our high redshift clusters, but if 
we assume that 20% of the light that would be on the 
red sequence at z = ends up in the BCG-I-ICL this 
would move jcrs,pred for the fourth model even further 
into agreement with local values, implying that no new 
processes are needed to reconcile the red sequence lumi- 
n osities of cluster s at 0. 4 < z < 0.8 with those locally. 

IGonzalez et al.l ()2007t ) find that the fraction of light in 




Fig. 11. — A comparison of j for galaxies in EDisCS clusters 
as computed by integrating the measured LF and by integrating 
the Schechter function fits to each cluster. The left-hand column 
is for red galaxies and the right-hand column is for blue galaxies. 
The open squares show the measured values with error bars on 
each quantity for the integral of LFzp and the open triangles (only 
shown in the top panels - the g-band) are for the integral of LFss- 
The vertical error bars are usually smaller than the points. The 
solid triangles on the bottom of each panel show the clusters with 
poor Schechter function fits. The solid diagonal line shows the one- 
to-one relation. The dotted and dashed lines show the least squares 
fit to the data for LFzp and LFss respectively. The measured LFs 
seem to have converged for all the red LFs and for the blue LFzp. 
We are missing some light for the faintest clusters when integrating 
LFss- 

the BCG-I-I CL decreases with i ncreasing cluster velocity 
dispersion. iLin fc MohH (|2004( l measure a BCG magni- 
tude that may include some ICL and they find that the 
luminosity fraction also decreases with increasing cluster 
mass. In both cases the trends have a large scatter at ve- 
locity dispersions corresponding to our clusters and it is 
possible that a trend is present in our data but masked 
by the large scatter within our own sample. If such a 
trend exists it is possible that the BCG and ICL build 
up at different rates and with a different response to the 
accretion history of the cluster. Hopefully more progress 
will be made with future high resolution simulations of 
clusters in a cosmological context, and direct measure- 
ments of the high redshift ICL. 

7. SUMMARY AND CONCLUSIONS 

In this paper we have measured the rest-frame opti- 
cal LFs for 16 clusters at 0.4 < z < 0.8 that are drawn 
from the ESQ Distant Cluster Survey (EDisCS). These 
clusters have a range in velocity dispersions and, in con- 
trast to massive x-ray selected high redshift clusters, are 
the progenitors of "typical" clusters in the local Universe 
(|Milvang- Jensen et al.ll2008l ). 

We determined membership for our clusters using a 
photometric-rcdshift-based technique and one based on 
statistical background subtraction. From a detailed com- 
parison of these two methods we concluded that the LF 
could only be robustly determined for red sequence galax- 
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Fig. 12. — A comparison of M200 

in EDisCS clusters to 

that from clusters in the SDSS, in the rest-frame g, r, and i-bands. 
The large solid squares show the EDisCS clusters, where the hori- 
zontal error bars show the 68% confidence intervals on M200 that 
stem from the uncertainty in the velocity dispersions. The small 
dots show the values for the SDSS clusters. The large circles are 
the geometric mean of the individual SDSS j'crs values in different 
mass bins. The solid vertical error bars on the SDSS points show 
the geometric standard deviation in each mass bin. The diagonal 
black line shows the slope that galaxies will lie on if they have con- 
stant jcrs/M200- The SDSS clusters deviate significantly from this 
relation indicating a real trend of jcrs/M200 with M200- 

ies and that the two methods resulted in very different 
M* and a values for blue galaxies. We therefore focus on 
the LFs of red sequence galaxies. We computed individ- 
ual LFs for our clusters and composite LFs for the whole 
sample as well as for sample split into subsets by redshift 
at z = 0.6 and velocity dispersion at a = 600 km/s. For 
the individual and composite LFs we fit Schechter func- 
tions, where we fixed the faint-end slope to the value 
determined from the fit to the composite of all EDisCS 
clusters in two different redshift bins. 

As a low redshift comparison sample we used a cluster 
catalog drawn from the SDSS and calculate the compos- 
ite LF and its Schechter function fit as for the EDisCS 
clusters. When splitting the SDSS sample into bins of ve- 
locity dispersion we take into account the average mass 
growth in clusters as expected from numerical simula- 
tions. In this way we can compare clusters at high red- 
shift to (representatives of) their likely descendants at 
low redshift, something that has not been possible with 
previous LF studies that either concentrated on very 
massive high redshift clusters - whose descendants would 
be largely absent from local volumes - or that had no 
velocity dispersion information for the clusters. 

We measure significant evolution in the LF of cluster 
red sequence galaxies at z < 0.8. In detail the LFs show 
evolution in the bright end consistent with passive evolu- 
tion but show a dramatic increase in the number of faint 
galaxies relative to bright ones toward lower redshifts, 
both within our own survey and when compared to the 
SDSS cluster sample. As a simple characterization of this 
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Fig. 13. — A comparison of jcrs/M200 vs M200 for red sequence 
galaxies in EDisCS clusters to that from clusters in the SDSS, in the 
rest-frame g, r, and i-bands. {left column) The large solid squares 
show the EDisCS clusters, where the vertical and horizontal error 
bars show the 68% confidence intervals on jcrs/M200 and the errors 
in M200 that stem from the uncertainty in the velocity dispersions, 
respectively. The small dots show the values for the SDSS clusters. 
The large circles are the geometric mean of the individual SDSS 
jcrs/M2oo values in different mass bins. The solid vertical error 
bars show the geometric standard deviation. Note that the errors 
m Jcrs /M200 and M200 are correlated. The arrow in the upper left 
panel demonstrates the change in jcrs/M2oo that results from a 
factor of four change in M2oo- {right column) The circles and error 
bars are the same as in the left column. The large open pentagons 
show the geometric mean of the EDisCS clusters in three mass 
bins chosen to contain roughly equal numbers of clusters, with the 
vertical error bars showing the geometric standard deviation for 
each value. The horizontal position of each pentagon is determined 
by the geometric mean of the masses for the EDisCS clusters in that 
mass bin. The arrows show the expected evolution from z^ciust to 
^ = for a model in which the cluster mass increases by accretion 
but where no new galaxies are added to the red sequence and where 
those that are already present at the epoch of observation fade 
passively as SSPs with Zf = 2. 



evolutio n we measure the ratio of luminous to faint galax- 
ies as in iDe Lucia et all ()2007t ) and find similar results. 
We also measure the build-up of the red sequence as a 
more detailed function of magnitude and find tentative 
evidence for an evolving magnitude threshold brighter 
than which the LF is in place w.r.t. the local LF. It is 
not clear if this evolving magnitude threshold is in any 
way related to the evolving mass threshold seen in field 
samples, above which star formation is truncated (e.g. 
Bundy et al. 2006), or if it corresponds to a different 
cluster related quenching mechanism. Indeed, the late 
build-up of the faint red sequence in clusters may also 
be related to the increase in the SO fraction seen toward 
lower redshift. 

We perform Schechter function fits to our LFs and find 
significant evolution in a but no evolution in M*, despite 
finding that the measured LFs at the bright end are con- 
sistent with passive evolution. This highlights the com- 
plications of using M* as a measure of the evolution in 
the luminosity of the galaxy population as a whole when 
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Fig. 14. — The velocity dispersions for our 16 clusters as a 
function of redshift. A Spearman Rank test gives a 4% probability 
that sigma is correlated with redshift. 
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Fig. 15. — The fading as a function of time for simple stellar 
populations in the g, r, and i-bands. The zeropoint of the y-axis 
is arbitrary. The thick solid line and thick dashed lines show the 
evolution in magnitude for SSP models with Z=Zq and 2.5Zq 
respectively. The thin dotted line is the linear fit to the solar 
metallicity model. The 2.5Zq model is brighter in the mean by 
0.2-0.4 magnitudes at these ages and we have subtracted out the 
difference to highlight the similar slopes. At (t) > 1.4 Gyr the 
fading in magnitude is well approximated by a linear relation in 
magnitude and log(i). 
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Fig. 16. — The ratio of the predicted red sequence g-band 
luminosity, jcrs,predj of EDisCS clusters at 2 = to the measured 
jcis from SDSS clusters, as a function of cluster mass. We show 
the predicted values for 4 toy models described in the text. The 
squares show a model in which only mass is added to the clusters 
but in which no new galaxies are added to the red sequence and 
those that exist at the epoch of observation evolve passively. The 
triangles show what happens when galaxies are added onto the 
red sequence with the same Mtot/L and its evolution as cluster 
galaxies. The pentagons show a model in which no additional 
red galaxies are accreted into the cluster but in which the blue 
galaxies in the clusters at z = 0.6 are assumed to have constant 
SFRs prior to the epoch of observation but subsequently have their 
SFR truncated and evolve passively thereafter. The circles show a 
model in which both the blue galaxies in the clusters are added to 
the red sequence and in which red passively evolving galaxies are 
accreted from the field. In all models the error bars account for 
the dispersion in jcrs values in the EDisCS data. In the third and 
fourth models the error bars also account for the different ji,g for 
different membership techniques and for the range of possible SFH. 
In each mass bin the points have been offset in mass for clarity. 

a is also simultaneously undergoing strong evolution. In 
our case it must be that luminosities of the whole galaxy 
population are not evolving in lock step. 

We split our sample into two bins of velocity disper- 
sion and find only small differences in the detailed LFs, 
although we find the same result as DL07 that the ratio 
of luminous to faint galaxies is higher, and evolves more 
quickly, in clusters of higher velocity dispersion. We find 
indistinguishable luminous-to-faint ratios for SDSS clus- 
ters of different velo cit y di spersion, similar to that found 
bv fde Propris et all (|200l . 

We looked for radial trends by examining the EDisCS 
LFs computed at r < 0.75Mpc and r < 0.5Mpc and find 
no difference in either M* or a for red sequence galaxies. 
This comparison, however, is uncertain as the two radial 
bins are highly correlated. 

Using the field L F of red sequence galaxies measured 
from the NDWFS ()Brown et al.l [2007b and MUSYC we 
compared our cluster LFs to the coeval field LF for sim- 
ilarly selected galaxies. At 0.6 < z < 0.8 the field 
has more faint galaxies relative to bright ones than the 
clusters but at 0.4 < z < 0.6 this has reversed, with 
the clusters having more faint galaxies than the field. 
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This epoch is similar to that in which the buildup of 
the SO population in clusters starts to become significant 
(|Desai et al.ll2007l ). Combined with the more rapid evo- 
lution of -/ViumZ-^faint for high vclocity dispersion clus- 
ters, the different rates of evolution in the LFs imply 
that dense environments are more efficient than the field 
at adding galaxies to the red sequence at z < 1. These 
trends in the ratio of luminous to faint galaxies are re- 
flected in the Schechter function fits. At both redshifts 
the EDisCS LF has a more negative a than the field but 
a slightly fainter M*. While the former agrees with local 
cluster-field comparisons from 2dF, the latter disagrees. 
Discovering the cause of this discrepancy in the relative 
M* values will require larger samples of clusters at in- 
termediate redshift to increase the signal-to-noise of the 
LF at high luminosities. 

To constrain different mechanisms for building up the 
red sequence galaxy population we measure the total red 
sequence light in the EDisCS and quantify its evolution 
w.r.t. clusters from SDSS. Clusters at high redshift are 
overluminous compared to their likely local descendants. 
Once passive fading is accounted for it appears that the 
clusters are a factor of 1 — 3 underluminous compared to 
the local clusters that are their likely descendants. Since 
light traces stellar mass on the red sequence, this implies 
that the mass on the red sequence in clusters must grow 
by a factor of 1 — 3 at z < 0.8, with most of the growth 
occurring at the faint end of the LF where we directly 
witness strong evolution. This is a similar amount of 
growth in the red sequence as inferred from studies of the 
field LF (Bell ct al. 2004; Fabcr ct al. 2007; Brown^eLalJ 
l2007t ) and indicates that the additional transformation of 
blue galaxies to red galaxies in clusters may be modest. 
However, due to the significant uncertainties we cannot 
determine if the total amount of mass added depends 
on cluster velocity dispersion. Evidence for environmen- 
tal effects come predominantly in the dependence of the 
shape evolution of the LF on cluster velocity dispersion 
and in differences between the cluster and field. 

To explore what physical mechanisms may be driving 
the assembly of the red sequence we explore a set of sim- 
ple toy models that incorporate many of the processes 
that should add galaxies to the red sequence in clusters 
in a conservative manner. Accounting for all necessary 
processes we find that these models overpredict the light 
in local clusters. The model predictions can be recon- 
ciled with the data by a combination of three previously 
known processes: blue galaxies that have fallen in since 
z < 0.6 but are not on the red sequence at 2: = 0, attenu- 
ation of light on the red sequence by dust extinction, and 
the transfer of stars from galaxies to the diffuse cluster 
light via tidal stripping. 

The results presented in this paper were only made 
possible with a large sample of clusters that span a range 
of redshift and velocity dispersion, that have accurately 
measured velocity dispersions, and that have deep multi- 
wavelength photometry over a significant fraction of the 
virial radius. Our analysis was limited in large part by 
our reliance on photometric methods to isolate cluster 
members, and by the limited number of clusters in our 
sample. To improve upon the analysis several ingredi- 
ents are needed. First, we need more clusters over a large 
range of mass, to confirm the mass dependent assembly 
of the red sequence. Second, we need to move to larger 



radii so that we can probe out past i?2oo f^nd thereby in- 
clude many more galaxies in our sample. Third, we need 
increased spectroscopy of blue cluster members to allow 
a robust luminosity function determination for galaxies 
of all colors. This last element would also be assisted if 
we had wide fields as we would then be able to do statis- 
tical background subtraction using "field" samples with 
identical photometry and roughly co-spatial on the sky, 
thus bypassing many of the problems inherent in using 
external fields for background estimation. 
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Zeropoints of fits to red sequence 
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Note. — Zeropoints of V — I vs. 
/tot color magnitude relation are calcu- 
lated for spectroscopically defined non- 
starforming galaxies. These are defined 
where /to t= 0, which diff ers from the def- 
inition of lDe Lucia et al.l l[2007i ). 
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SDSS LF Schechter Function Parameters 
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< 


-22.0 


87.44 ± 15.39 


42.04 ± 9.32 


23.16 ±6.73 


18.15 ±5.93 


42.87 ± 11.48 


23.65 ± 7.69 


17.65 ±7.16 


-22 





< 


Mi 


< 


-21.5 


108.61 ± 16.03 


64.97 ± 10.42 


54.26 ±9.53 


15.84 ± 5.39 


41.78 ± 11.44 


23.35 ±7.12 


17.12 ± 7.36 


-21 


5 


< 


Mi 


< 


-21.0 


119.74 ± 17.39 


67.00 ± 10.95 


43.61 ± 8.72 


24.50 ± 6.54 


52.02 ± 12.84 


31.93 ± 8.93 


19.85 ± 7.98 


-21 





< 


Mi 


< 


-20.5 


97.43 ± 15.34 


68.10 ± 11.41 


31.05 ± 7.82 


33.68 ± 7.46 


29.81 ± 9.92 


17.13 ±6.08 


11.93 ± 6.41 


-20 


5 


< 


M, 


< 


-20.0 


88.31 ± 13.66 


72.39 ± 11.42 


39.31 ± 8.33 


31.64 ±7.23 


17.87 ±7.56 


17.38 ±7.66 


3.10 ± 3.45 


-20 





< 


Mi 


< 


-19.5 


78.96 ± 12.73 


68.22 ± 11.14 


39.00 ±8.16 


28.52 ±7.03 


12.40 ± 6.29 


8.59 ± 4.32 


4.01 ± 3.86 


-19 


5 


< 


Mi 


< 


-19.0 


51.95 ± 11.05 


42.79 ±9.16 


26.64 ± 7.08 


16.45 ±5.59 


10.27 ±6.18 


8.44 ± 5.24 


2.59 ± 3.40 


-19 





< 


Mi 


< 


-18.5 




51.44 ± 10.03 


32.30 ± 7.75 


19.60 ±6.12 






4.01 ± 3.86 


-18 


5 


< 


Mi 


< 


-18.0 
















-18 





< 


Mi 


< 


-17.5 

















CD 
P 

o 
O 



o 

A 
p 
bo 



Note. — Magnitudes are given in units of M — 5 log hro- Composite LFs are only given for magnitudes brighter than which all 
clusters in each redshift and velocity bin are complete. 



Rudnick et al. 



TABLE 4 

SCHECHTER FUNCTION PARAMETERS FOR EDisCS COMPOSITE LFS 



redshift 




M* 

M — 5 log hfo 


"9 


M - 5 log hro 




M* 

M — 5 log hro 




0.4 < z < 0.8 
0.4 < z < 0.6 
0.4 < z < 0.6 
0.4 < z < 0.6 
0.6 < z < 0.8 
0.6 < z < 0.8 
0.6 < z < 0.8 


all clusters 
all clusters 

> 600 km/s 

< 600 km/s 
all clusters 

> 600 km/s 

< 600 km/s 


-20.92t°?i 

-20.76l°;f, 

-20.73iO:lf, 

-20.80+014 

-20.79l°;« 

-20.77t°-l^ 
-20.81±0-ie 


b 
b 

a 
a 


-21 

Z1.01_o 14 

-21.48l°i« 

-21.38l«:i2 

-2i.5oi;;:fi 

— 91 41 +0-88 
"'-'-•^-'--().20 

-21 'W+O-IS 
^-l--'^' -0.14 

-21.46l«i2 


-n s6+°-i'' 
^■•^''-o.os 

b 
b 

o.o8+°:t^ 

a 

a 


-21.80lo.?? 
-21.83t°-28 

-21.8ltrJ 

01 7'c:+0.22 

-21.64l°j8 

-21.67t°-}| 

-2i.66iS!:?l 


'J-^^-o.io 

-U.bS^O 06 

b 

b 

1 7+0-42 
-0.20 

a 

a 



Uses distribution of a determined from fits to all clusters at 0.4 < z < 0.6.'' Uses distribution of a determined from fits 
to all clusters at 0.6 < z < 0.8. 



APPENDIX 

REST-FRAME OPTICAL LUMINOSITY FUNCTIONS OF EDISCS CLUSTERS 



TABLE 5 

Rest-frame ^-band LPs for EDisCS clusters at OA < z < 0.6 



Mg 

M - 5 log /i7o 



*cll018. 8-1211 *cll037.9-1243 ^cll059. 2-1253 *clll38. 2-1133 *cll202. 7-1224 *cll232.6- 



*oll301. 7-1139 *cll3B3. 0-1137 *cll411. 1-1148 ^cll420. 3-1236 



-24.5 < Mg < 
-24.0 < Mg < 
< 
< 
< 



-23.5 < Mg 
-23.0 < Mg 
-22.5 < Mn 



-22.0 < Mg < 

-21.5 < Mg < 
-21.0 < Mg < 

-20.5 <Mg < 

-20.0 <Mg < 
-19.5 < Mg < 

-19.0 < Mj < 

-18.5 < Mg < 
-18.0 <Mg < 
-17.5 < Mg < - 

-17.0 <Ma < 



24.0 
23.5 
23.0 
22.5 
22.0 
21.5 
21.0 
20.5 
20.0 
19.5 
19.0 
18.5 
18.0 
17.5 
17.0 
16.5 



< 2.00 

< 2.00 

< 2.00 

< 2.00 

< 2.00 

1+2.41 

-^-0.83 
1+2.41 
-^-0.83 

q + 3.00 

•^-1.63 
4+3.24 
^-1.91 
14+4.87 



-3.69 
c:+3.4S 
°-2.15 
7+3.83 
'-2.58 
c:+3.45 
°-2.15 
4+3.24 
^-1.91 
1+2.41 
^-0.83 



< 2.00 

< 2.00 

< 2.00 

< 2.00 

< 2.00 

1+2.41 
^-0.83 
c+3.45 
^-2.15 
O+3.00 
•^-1.63 
4+3.24 
^-1.91 
9+2.73 
^-1.29 
c+3.45 
^-2.15 
4+3.24 

O+3.00 
■'-1. 63 

< 2.00 

Y+3.83 
'-2.58 



< 2.00 

< 2.00 

< 2.00 

< 2.00 

1+2.41 
^-0.83 

c+3.45 
^^-2.15 
7+3.83 
'-2.58 
1Q+4.32 

q+4.16 
^-2.94 
Q+4.16 
="-2.94 

c+3.65 
"-2.37 
O+4.00 
"-2.76 
O+3.00 
■'-1. 63 
c:+3.45 
^-2.15 
9+2.73 
^-1.29 



< 2.00 

< 2.00 

< 2.00 

< 2.00 

< 2.00 

1+2.41 
^-0.83 



O+3.00 
■'-1. 63 
1Q+4.32 

q+4.16 
^-2.94 
fi+3.66 
"-2.37 

O+4.00 
»-2.76 
(5+3.65 
°-2.37 
4+3.24 
^-1.91 
c+3.45 
°-2.15 
9+2.73 
^-1.29 
9+2.73 
^-1.29 



< 2.00 

< 2.00 

< 2.00 

< 2.00 

1+2.41 
^ -0.83 



1+2.41 

-'--o.ss 
< 2.00 

1+2.41 
'--0.83 
c+3.45 
^-2.15 
fi+3.66 
"-2.37 
O+4.00 
*'-2.76 
9+2.73 
^-1.29 
1+2.41 
'-0.83 
9+2.73 
^-1.29 
1+2.41 

-'■-o.ss 



< 2.00 

< 2.00 

< 2.00 

1+2.41 
^-0.83 
Q+3.00 
•^-1.63 

4+3.24 
^-1.91 
^3+4.74 



16-f 



''-3.55 
,+5.12 
-3.95 
1 9+4.61 
^^-3.41 
19+4.61 
^"^-3.41 

^g + 5.47 

^=^-4.32 
19+4.61 
-^^-3.41 
19+4.61 
-^^-3.41 
94+6.00 
^^-4.86 



< 2.00 

< 2.00 

< 2.00 

< 2.00 

< 2.00 

1+2.41 
-'--0.83 
7+3.83 
' -2.58 
Q+4.16 
"-2.94 
Q+3.00 
•^-1.63 
fi+3.65 
"-2.37 
3+3.00 



1.63 
+3.65 
2.37 
+3.83 
-2.58 
+4.46 
3.26 
+4.32 
3.10 



< 2.00 

< 2.00 

< 2.00 

< 2.00 

9+2.73 
^-1.29 
1+2.41 
^-0.83 
4+3.24 
^-1.91 
f;+3.65 
"-2.37 
4+3.24 
^-1.91 
g+3.65 



-2.37 
4+3.24 
^-1.91 
4+3.24 
^-1.91 
fi+3.65 
"-2.37 
4+3.24 
^-1.91 



< 2.00 

< 2.00 

< 2.00 

< 2.00 

1+2.41 
^-0.83 
9+2.73 
^-1.29 
1+2.41 



6 



■--0.83 
+3.65 
-2.37 
ic+5.00 
^"^-3.83 
fi+3.65 
"-2.37 
Q+4.16 
"-2.94 
O+4.00 
°~2.76 
1+2.41 
-'-0.83 
Q+3.00 
•^-1.63 



< 2.00 

< 2.00 

< 2.00 

< 2.00 

9+2.73 
^-1.29 
1+2.41 
^-0.83 



Q+3.00 
■'-1.63 
7+3.83 
'-2.58 
Q+3.00 
•^-1.63 
r;+3.46 
^-2.15 
fi+3.65 
"-2.37 
Q+3.00 
■'-1. 63 
Q+3.00 
•'-1. 63 
4+3.24 
^-1.91 
1+2.41 
-'■-0.83 



Note. — $ gives the number of galaxies per magnitude bin for ea^h cluster and does not depend on h-ro- 



TABLE 6 

Rest-frame ^-band LPs for EDisCS clusters at 0.6 < z < 0.8 



*cll040.7-llBB *cll064. 4-1146 *cll054. 7-1245 *oll216. 8-1201 *cll227.9-1138 *cll3B4.2-1230 

M - 5 log ?170 



-24.5 


< 


Mg 


< 


-24.0 


< 2.00 


< 2.00 


< 2.00 


< 2.00 


< 2.00 


< 2.00 


-24.0 


< 


Mg 


< 


-23.5 


< 2.00 


< 2.00 


< 2.00 


< 2.00 


< 2.00 


< 2.00 


-23.5 


< 


Mg 


< 


-23.0 


< 2.00 


< 2.00 


< 2.00 


< 2.00 


< 2.00 


< 2.00 


— zo.U 


< 


Mg 




— zz.o 


< Z.UU 


9 nn 


1+2.41 
-^-0.83 


Q+3.00 
•^-1.63 


< z.UU 


^ 9 An 
< z.UU 


—22.5 


< 


Mg 


< 


—22.0 


9+2.73 
^-1.29 


O+3.00 
•^-1.63 


1 +2.41 
-^-0.83 


4+3-24 
^-1.91 


< 2.00 


1+2.41 
^-0.83 


—22.0 


< 


Mg 


< 


—21.5 


n + 3.00 

■"'-1.63 


„-f-3,00 
■'-1. 63 


„-f-3.00 
■^-1.63 


Q+4 74 


< 2.00 


Q+3.00 
•^-1.63 


-21.5 


< 


Mg 


< 


-21.0 


(5+3.65 
"-2.37 


1 1 +4.46 
^^-3.26 


O+4.00 
"-2.76 


1Q+4.32 
^"-3.10 


9+2.73 
^-1.29 


9+2.73 
^-1.29 


-21.0 


< 


Mg 


< 


-20.5 


O+3.00 
■'-1. 63 


O+4.00 
"-2.76 


1Q+4.32 
J^'J-3.10 


14+4.87 
^^-3.69 


Q+3.00 
■^-1.63 


Q+3.00 
•^-1.63 


-20.5 


< 


Mg 


< 


-20.0 


Q+3.00 
•^-1.63 


11+4.46 
^^-3.26 


q+4.16 
^-2.94 


14+4-87 
^^-3.69 


4+3.24 
^-1.91 


4+3-24 
^-1.91 


-20.0 


< 


Mg 


< 


-19.5 


9+2.73 
^-1.29 


1Q+4.32 
-^"-3.10 


Y+3.83 
'-2.58 


O+4-00 
°-2.76 


9+2.73 
^-1.29 


Q+3.00 
•^-1.63 


-19.5 


< 


Mg 


< 


-19.0 


< 2.00 


9+2.73 
^-1.29 


7+3.83 
'-2.58 


9+4.16 
"-2.94 


< 2.00 


9+2.73 
^-1.29 


-19.0 


< 


Mg 


< 


-18.5 


1+2.41 
-^-0.83 


y+3.83 

'-2.58 


O+4.00 
"-2.76 
tr+3.45 
"-2. 15 


-,Q+4.32 
^"-3.10 
f^+3.65 
"-2.37 


< 2.00 


1+2.41 
^-0.83 


-18.5 


< 


Mg 


< 


-18.0 


< 2.00 


0+3.00 
■'-1.63 


< 2.00 


< 2.00 


-18.0 


< 


Mg 


< 


-17.5 


9+2.73 
^-1.29 


Q+3.00 
•^-1.63 


9+4.16 
^-2.94 


a+4.00 
°-2.76 


< 2.00 


Q+3.00 
•^-1.63 


-17.5 


< 


Mg 


< 


-17.0 










1+2.41 
^-0.83 




-17.0 


< 


Mg 


< 


-16.5 















Note. — $ gives the number of galaxies per magnitude bin for each cluster and does not depend on hjo ■ The LF for 
CL1227.9-1138 has been computed over < 50% of the full cluster area and so must be re-normalized by the full area. 



TABLE 7 

Rest-frame r-BAND LFs for EDisCS clusters at 0.4 < z < 0.6 



M ■ 



Mr 

5 log h7o 



*cll018. 8-1211 *cll037.9-1243 *cll059.2-1263 *clll38. 2-1133 ^cll202. 7-1224 *cll232.5- 



*cll301. 7-1139 *cll3B3. 0-1137 *cll411. 1-1148 *cll420.3-1236 



-24, 


.5 


< 


Mr 


< 


-24, 


,0 


— 24, 


.0 


< 


Mr 




—23, 


.5 


— Zo. 


c 
.O 


< 


iWy 




— Zo. 


.U 


— Zo. 


.u 


< 


iVlr 




— zz. 


.0 


—22, 


.5 


< 


Mr 




—22, 


.0 


-22, 


.0 


< 


Mr 


< 


-21, 


,5 


-21, 


.5 


< 


Mr 


< 


-21, 


,0 


-21, 


.0 


< 


Mr 


< 


-20. 


.5 


-20, 


.5 


< 


Mr 


< 


-20. 


.0 


-20, 


.0 


< 


Mr 


< 


-19, 


,5 


-19, 


.5 


< 


Mr 


< 


-19, 


,0 


-19, 


.0 


< 


Mr 


< 


-18, 


.5 


-18, 


.5 


< 


Mr 


< 


-18. 


.0 


-18, 


.0 


< 


Mr 


< 


-17. 


.5 


-17, 


.5 


< 


Mr 


< 


-17, 


,0 


-17 


,0 




Mr 


< 


-16 





< 2.00 

< 2.00 

< 2.00 

< 2.00 

r, + 2.73 

^-1.29 

q + 3.00 

"^-1.63 
2+2.73 



1.29 
+3.83 
2.58 
+4.46 
3.26 
+3.00 
1.63 
+4.16 
2.94 
+3.24 
1.91 
+3.00 
1.63 



r 
11 

3 

9 

4; 

3 

< 2.00 



< 2.00 

< 2.00 

< 2.00 

< 2.00 

r, + 2.73 

^-1.29 
4+3-24 
^-1.91 
4+3-24 
^-1.91 

q + 3.00 

•^-1.63 

q+3.00 

•^-1.63 
4+3.24 
^-1.91 
c+3.45 
0-2.16 

r, + 2.73 

^-1.29 
4+3.24 

r, + 2.73 

^-1.29 



< 2.00 

< 2.00 
1+2.41 

-Q-83 

J+3, 

10 
T 
10 
T 

6 



+2.73 
1.29 
" 24 
91 
+4.32 
3.10 
+3.83 
2.58 
+4.32 
3.10 
+3.83 
2.58 
+3.65 
2.37 
+4.00 
-2.76 
+3.00 
-1.63 
+3.45 
-2.15 
+2.73 
-1.29 



< 2.00 


< 2.00 


< 2.00 


< 2.00 


< 2.00 


< 2.00 


< 2.00 


< 2.00 


<. z.UU 


^ o nn 


-1+2.41 
-^-0.83 


^ o nn 
< z.uu 


-1+2.41 
-^-0.83 


-1+2.41 
^-0.83 


4+3.24 
4-1.91 


1+2.41 
-^-0.83 


< 2.00 


-1+2.41 
^-0.83 


c+3.45 
0-2.15 


1+2.41 
-^-0.83 


Q+3.00 
•^-1.63 


< 2.00 


1 q+4.74 


O+4.00 
°-2.76 


iq+4.74 


r, + 2.73 

^-1.29 


-17+5.24 
^'-4.08 


Q+4.16 
''-2.94 


O+4.00 


c+3.45 


,„+4.32 
J^^-3.10 


9+2.73 


°-2.76 


0-2.15 


^-1.29 


ir+3.46 
0-2.15 


(;+3.65 


19+4.61 
-■-^-3.41 


7+3.83 
'-2.58 


°-2.37 


O+4.00 


O+4.00 


lQ+5-47 
J^^-4.32 


9+2.73 
^-1.29 


°-2.76 


°-2.76 


ir+3.46 
0-2.15 


9+2.73 
^-1.29 


-1-1+4.46 
^^-3.26 


4+3.24 
^-1.91 


4+3.24 
^-1.91 


< 2.00 


17+5.24 
^'-4.08 


7+3.83 
'-2.58 


Q+3.00 
•^-1.63 


Q+2.73 
^-1.29 


iq+6-47 
-■■^-4.32 


19+4.61 
^^-3.41 


ir+3.46 
0-2.15 


< 2.00 




,(,+4.32 
^^^-3.10 


Q+2.73 

^-1.29 


-1+2.41 











< 2.00 

< 2.00 

-1+2.41 
^-0.83 
1+2.41 
^-0.83 
1+2.41 
^-0.83 
7+3.83 
'-2.68 
2+2.73 



1.: . 
c:+3.45 
0-2.16 
fi+3.65 
"-2.37 
4+3.24 
^-1.91 
4+3.24 



91 
+3.00 
-1.63 
+3.00 
-1.63 



< 2.00 

< 2.00 

< 2.00 

1 

2 

5; 

5 

13 
5 

10 
6 
2 



+2.41 
0.83 
+2.73 
1.29 
+3.45 
2.16 
+3.45 
2.15 
-4.74 
— 3.55 
+3.45 
2.15 
-4.32 
-3.10 
+3.65 
2.37 
+2.73 
1.29 



,+2.73 
^-1.29 
1+2.41 
■^-0.83 



< 2.00 

< 2.00 

< 2.00 

9+2.73 
^-1.29 

q + 3.00 

"^-1.63 
4+3.24 
^-1.91 
4+3.24 
^-1.91 

q+3.00 

•^-1.63 
O+4.00 
°-2.76 
4+3.24 
*-1.91 
9+2.73 
^-1.29 

q + 3.00 

■^-1.63 
4+3.24 
*-1.91 
1+2.41 
^-0.83 



Note. — $ gives the number of galaxies per magnitude bin for each cluster and does not depend on h-ro- 



to 

CO 



CO 

TABLE 8 

Rest-frame t-band LFs for EDisCS clusters at 0.6 < 2; < 0.8 



Mr *cll040. 7-1166 *cllOB4.4-1146 *cllOB4.7-1246 *cll216. 8-1201 *cll227.9-1138 *cll364. 2-1230 

M — 5 log ^70 

< 2.00 < 2.00 < 2.00 < 2.00 < 2.00 < 2.00 

< 2.00 < 2.00 < 2.00 lto83 < 2-00 < 2.00 
<2.00 itltl itltl 2till <2-00 <2.00 

.+3.24 r,+2.73 1+2.41 o+4.00 ^ 9 nO 1+2.41 

^-1.91 ^-1.29 -^-0.83 °-2.76 ^ "^-0.83 

O+3.00 O+4.00 4+3.24 -,9+4.61 -,+2.41 ^+3.24 

■^-1.63 "-2.76 ^-1.91 -^^-3.41 -^-0.83 *-1.91 

C-+3.45 Q+4.16 Q+4.16 -|r|+4.32 r,+2.73 9+2.73 

''-2.15 ^-2.94 ^-2.94 -"-^-3.10 ^-1.29 ^-1.29 

9+2.73 Q+4.16 ir,+4.32 ,9+4.61 o+3.00 4+8. 24 

^-1.29 ^-2.94 ^"-3.10 ^^-3.41 "^-1.63 ^-1.91 

4+3.24 ,f,+4.32 0+4.00 1J-+5.00 /I+3.24 9+2. 73 

*-1.91 ^"-3.10 °-2.76 ^''-3.83 ^-1.91 ^-1.29 

-,+2.41 0+4.00 7+3.83 7+3.83 -,+2.41 9+2.73 

-"--0.83 °-2.76 '-2. ,■38 ' -2.,S8 -'--0.83 ^-1.29 

< 2.00 2t-g 6t^:« 8t^:?» < 2.00 3l?:« 

-1+2.41 4+3.24 c:+3.45 7+3.83 ^ 9 nn -1+2.41 PJ 

-^--0.83 ^-1.91 ^-2.15 -2.58 ^ z.uu -^-0.83 5 

<2.00 4tl l* ^ii'gi ^ii'gi < 2-00 < 2.00 ^ 

itltl <2.00 <2.00 ... ^ 

p 
o 

Q 

Note. — $ gives the number of galaxies per magnitude bin for ea<;h cluster and does not depend on hro ■ The LF for 
CL1227.9-1138 has been computed over < 50% of the full cluster area and so must be re-normalized by the full area. S- 
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TABLE 9 

Rest-frame i-band LFs for EDisCS clusters at 0.4 < 2 < 0.6 



M- 



Mi 

5 log h70 



*cll018. 8-1211 *cll037.9-1243 *cll069. 2-1253 *clll38. 2-1133 *cll202.7-1224 *cll232.6-1260 *cll301. 7-1139 *cll3B3. 0-1137 *cll411. 1-1148 *cll420. 3-1236 
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< 2.00 

< 2.00 

< 2.00 
1 
1 
3 
5 
10 
T 
7 
6 



+2.41 
-0.83 
+2.41 
-0.83 
+3.00 
1.63 
+3.45 
2.16 
+4.32 
3.10 
+3.83 
2.58 
+3.83 
-2.58 
+3.65 
2.37 
4+3.24 
^-1.91 

< 2.00 



< 2.00 

< 2.00 

< 2.00 

1+2.41 
-0.83 
-+3.45 
-2.15 
:,+2.73 
-1.29 
,+3.00 
-1.63 
4+3.24 
^-1.91 
9+2.73 
^-1.29 
rr+3.45 
^-2.15 
9+2.73 
^-1.29 
Q+3.00 
•^-1.63 

< 2.00 



< 2.00 

< 2.00 

1+2.41 
-0.83 
.1+3.24 
-1.91 
a+4.00 
-2.76 
f3.83 
-2.58 
□+4.16 
-2.94 
D+4.00 
-2.76 
7+3.83 
-2.58 
s+4.00 
-2.76 
-+3.45 
-2.15 
-+3.45 
-2.15 
,+2.73 
-1.29 
1+2.41 
-^-0.83 
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< 2.00 
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^-1.29 
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7"f 
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ir+3.46 
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O+4.00 
°-2.76 
c+3.45 
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c+3.45 
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•^-1.63 
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*-1.91 
9+2.73 
^-1.29 



< 2.00 
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1+2.41 
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1+2.41 
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< 2.00 

4^ 

7l 
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'-2.58 
4+3.24 
^-1.91 
1+2.41 
^-0.83 
9+2.73 
^-1.29 

< 2.00 

1+2.41 

-'■-o.ss 



^+3.24 
-1.91 

7+3.83 
-2.58 

7+3.83 
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< 2.00 
1+2.41 
^-0.83 
q+3.00 
•^-1.63 
O+3.00 
■'-1.63 
+4.32 
3.10 
17+5.24 
' -4.08 
14+4-87 
^^-3.69 
+3.83 
-2.58 
+5.12 
-3.95 
+5.24 
-4.08 
+4.87 
3.69 



7Z 

17+ 
14+ 

1 1 +4:46 
^^-3.26 
+4.32 
3.10 



10 



< 2.00 

< 2.00 

< 2.00 

1+2.41 
-^-0.83 
c:+3.45 
^-2.15 
11+4.46 
^^-3.26 
4+3.24 
^-1.91 
4+3.24 
^-1.91 
n+3.00 
•^-1.63 
f^+3.65 
"-2.37 
9+2.73 
^-1.29 
19+4.61 
-^^-3.41 
7+3.83 



< 2.00 

< 2.00 

,+2.73 
-1.29 
1+2.41 

-0.83 
.,+3.24 
-1.91 
-+3.45 
-2.15 
4+3.24 
^-1.91 
-+3.45 
-2.15 
-+3.45 
-2.15 
5+3.00 
-1.63 
.,+3.24 
-1.91 
4+3.24 
^-1.91 



< 2.00 

< 2.00 

+2.41 
0.83 
9+2.73 
^-1.29 
+2.41 
0.83 
+3.65 
2.37 
+4.74 
3.55 
+3.65 
2.37 
+3.83 
2.58 
+4.16 
2.94 
+3.00 
1.63 
1+2.41 
^-0.83 
Q+3.00 
•^-1.63 



< 2.00 

< 2.00 

9+2.73 
^-1.29 
2+2.73 



-1.29 
Q+3.00 
■^-1.63 
c+3.45 
^-2.15 
Q+3.00 
•^-1.63 
Q+3.00 
•^-1.63 
O+4.00 
"-2.76 
4+3.24 
^-1.91 
9+2.73 
^-1.29 
4+3.24 
^-1.91 
9+2.73 
^-1.29 



Note. — $ gives the number of galaxies per magnitude bin for each cluster and does not depend on h-ro- 
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TABLE 10 

Rest-frame i-band LFs for EDisCS clusters at 0.6 < z < 0.8 



M ■ 



Mi 

5 log h70 



*cll040. 7-1166 *cllOB4.4-1146 *cllOB4.7-1246 *cll216. 8-1201 *cll227.9-1138 *cll364. 2-1230 
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9+2.73 
^-1.29 


Q+3.00 
•^-1.63 


9+2.73 
''-1.29 


O+3.00 
■^-1.63 


4+3.24 
^-1.91 


Q+3.00 
"^-1.63 


/^+3.65 
"-2.37 


-,(,+4.32 
-^^-3.10 


7+3.83 
'-2.58 


O+3.00 
•^-1.63 


1Q+4.32 


9+4.16 
^-2.94 


Q+3.00 
•^-1.63 


9+4.16 
^^-2.94 


9+4.16 
^^-2.94 


9+2.73 
''-1.29 


9+4.16 
="-2.94 


7+3.83 
'-2.58 


< 2.00 


9+2.73 
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fi+3.65 
"-2.37 


1+2.41 
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Q+3.00 
•^-1.63 
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< 2.00 

Q+3.00 
"^-1.63 
4+3.24 
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2.94 
^4.74 
— 3.55 
1C+5.00 
^"^-3.83 
Q+4.16 
="-2.94 
7+3.83 
'-2.58 
fi+3.65 
"-2.37 
fi+3.66 
"-2.37 



^ 9 on 
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9+2.73 
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4+3.24 
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9+2.73 
^-1.29 
Q+3.00 
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4+3.24 
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Note. — $ gives the number of galaxies per magnitude bin for each cluster and does not depend on hjo ■ The LF for 
CL1227.9-1138 has been computed over < 50% of the full cluster area and so must be re-normalized by the full area. 



